UNCLASSIFIED 


_ AD  NUMBER _ 

ADB016788 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  AUG  1976.  Other 
requests  shall  be  referred  to  Air  Force 
Armament  Lab.,  Eglin  AFB,  FL. 


_ AUTHORITY 

USADTC  ltr  10  Dec  1979 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 


AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200.20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  J°0N 
ITS  USE  AND  DISCLOSURE. 


DISTRIBUTION  STATEMENT  A 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


DDC  FILE  COPY _  ADB016788 


N 


AFATL-TR-76-90  V  y _ 

) 

DESIGN  CONSIDERATIONS 

FOR  AN  IMPULSE  CORRECTED,  BODY  FIXED, 

LASER  GUIDED,  FIN  STABILIZED  ROCKET 


BALLISTICS  BRANCH 

GUNS,  ROCKETS  AND  EXPLOSIVES  DIVISION 


AUGUST  1976 


y 


FINAL  REPORT:  JANUARY  1974  -  JULY  1976 


Distribution  limited  to  U.  S.  Government  agencies  only- 
tins  report  documents  test  and  evaluation;  distribution 
limitation  applied  August  1976  .  other  requests  for 
this  document  must  be  referred  to  the  Air  Force  Armament 
Laboratory  (DLDL),  hglin  Air  Force  Base,  Florida  32542. 


AIR  FORCE  ARMAMENT  LABORATORY 

AIR  FORCE  SYSTEMS  COMMAND  •  UNITED  STATES  AIR  FORCE 


EGLIN  AIR  FORCE  BASE,  FLORIDA 


UNCLASSIFIED 


SECURITY  CLASS'FICATION  OF  THIS  “AGE  ('Wh.ii  Bata  Hr, luted) 


REPORT  DOCUMENTATION  PAGE 


.  aEBiaaT.il  gi  in«» 


7  AUTHORS 

V 

Kenneth  K./Cobb 


AFATL-TR-76-9, 


n 


1/ 


2.  GOVT  ACCESSION  NO 


4.  TITLE  Canri  Subtitle) 


RESIGN  CONSIDERATIONS  FOR  AN , IMPULSE ^CORRECTED, 
JODY  FIXED,  LASER  GUIDED,  FIN  STABILIZED  ROCKET. 

^  !:  __  _  1l_  11.  _1_-H _ 


•A 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Guns,  Rockets  and  Explosives  Division  (DLDL) 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Armament  Laboratory 
Armament  Development  and  Test  Center 
Eglin  Air  Force  Base,  FL  32542 _ 


14.  MONITORING  AGENCY  name  »  ADORESSf/f  different  from  Controlling  Office) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


E Cl°’ FF T'S  CATALOG  NUMBER 


tv."  cr°rihT  A  PERIOD  COVERED 

, FINAL  'ye.f’t* 
lana—r  H74- 


PCArORMIRB  8RB. 


J  A  PERIOD  COVERED 

”  \ 
-Julg  tf76  a  \ 

IB.  RePORT  rBmber 


8.  CONTRACT  OR  GRANT  NUMBERfaj 


10  PPOGPAM  F.LEMENT,  PROJECT,  TASK 
APEA  A  WOF -  - 


6260 2  F 
25471)308 


ORjA-tTN.IT  NUMBERS 

OmM 


12.  REPORT  DATE 


Augil  «7(T7 


13.  NUMBER  OF  PAGES 
101 


13.  SECURITY  CLASS,  (of  Ihle  report) 

UNCLASSIFIED 


15a.  OECLASSl  FI  CATION/' DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Distribution  limited  to  U.S.  Government  agencies  only;  this  report  documents 
test  and  evaluation;  distribution  limitation  applied  August  1976.  Other  requests 
for  this  document  must  be  referred  to  the  Air  Force  Armament  Laboratory  (DLDL) , 
Eglin  Air  Force  Base,  Florida  32542. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetract  entered  In  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


Available  in  DDC 


19.  KEY  WORDS  (Continue  on  reverse  aide  it  necessary  and  identify  by  block  number) 


Body  Fixed  Guidance  System 
Impulse  Correction 


20.  ABSTRACT  (Continue  on  reverse  side  if  necessary  and  Identity  by  block  number) 


This  report  is  an  analysis  of  the  dynamics  and  accuracy  of  a  laser  guided  fin 
stabilized  rocket.  The  guidance  uses  a  quadrant  detector,  pursuit  guidance, 
and  is  a  body  fixed  system  (no  gyroscopic  or  relative  wind  reference).  Cor¬ 
rections  are  made  using  short  action  time  explosive  strips  located  forward  of 
the  center  of  mass. 

Several  guidance  consideration  which  affect  the  accuracy  or  effectiveness  of  — ’ 


-v" 

t-J 

/ 


DD 


FORM 
1  JAN  73 


1473  EDITION  OF  1  NOV  6S  IS  OBSOLETE 


UNCLASSIFIED 


40O 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (ItTien  Data  Entered)  ^ 


**-  C4 


SECUKITI  CLASSIFICATION 


JiflHiHmillgM 


(item  : 
cons  id' 
quency 
the  ro 
and  ae 
dead  z 


concluded)  This  system  are  analyzed  in  this  repoi 
at ions  are-  phasing  of  the  laser  pulse  frequency 
the  relation  between  the  trim  misalignment,  the  p 
rate;  the  maximum  expected  correction  as  a  funct 
dynamic  characteristics  of  the  rocket;  target  mot 

e  size. . 


Some  of  these 
the  pitch  fre- 
frequency,  and 
of  the  physical 
and  wind;  and  the 


m au* 


PREFACE 


tz\ «»“*»  B«„ch  of 

program  in  July  1973  The  knowloHo,.  'j  lhe  effort  commenced  on  this 
«U  continue  to  be  us J  to  supJo«  th^Ll  ^  ihiS."°r,!  has  bae"  used  and 
The  effort  has  been  carried  on  oerioH^.n  Tactical  Rocket  System, 

debted  to  Billy  Walker  and  Gus  Osterlow  fol  ^nce  1973<  1716  author  is  in- 

cir“^r 

assistance  in  computer  numerical  analysi ITsTj  IT^lr  p^X- 

TV*  i  ~  a.  _  . 


x - vgi  ajiuiixng . 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

Tlir 


FOR  THE  COMMANDER 


GERALD  P.  D'ARCV  ColinS?  ™ 

Chief,  Guns,  Rockets  and  Exp&lives  Division 


JSAF 


(The  reverse  of  this  page  is  blank) 


jliiiit  ii  ~  i  an  ■  mi 


I 


i  1 

I 


TABLE  OF  CONTENTS 


Section  Title  Page 

I  INTRODUCTION .  1 

II  GENERAL  GUIDANCE  CONSIDERATIONS  AND  SIMULATION  MODELS  ...  5 

1.  Basic  Guidance  System  .....  .  5 

2.  Nominal  Trajectory,  Baseline  Rocket  .  9 

3.  Linear  Pitching  and  Swerve  Motion  .  13 

4.  Six-Degree-of-Freedom  Guidance  Simulation .  20 

III  RELATION  BETWEEN  THE  LASER  PULSE  FREQUENCY  AND 

PITCH  FREQUENCY .  27 

IV  APPROXIMATE  MAXIMUM  EXPECTED  CORRECTION  .  36 

1.  Correction  Assuming  Impulse  Exactly  Toward  Target  ...  36 

2.  Effect  of  Four  Quadrant  Detector .  39 

3.  Maximum  Expected  Correction  for  Nominal 

Configurations .  42 

V  TRIM  MISALIGNMENT  EFFECT .  45 

VI  OTHER  GUIDANCE  CONSIDERATIONS  AND  ACCURACY  ERROR  SOURCES.  .  52 

1.  Parameters  Affecting  Accuracy .  52 

2.  Estimate  of  the  Accuracy  of  the  System .  56 

VII  DISTURBED  FLOW  EFFECT  AND  TEST .  60 

VIII  SUMMARY  AND  CONCLUSIONS .  67 

References  69 

Appendix 

A  ANGLE  OF  ATTACK  AND  SWERVE  FOR  PLANAR  MOTION,  ROLL 

RATE,  AND  ROLL  AVERAGED  THRUST .  71 

B  LINEARIZED  SOLUTION  FOR  THE  MOTION  OF  A  FIN- STABILIZED 

MISSILE  WITH  TRIM  MISALIGNMENT .  81 

C  FORCES  AND  MOMENTS .  85 

1 

D  EQUIVALENT  CEP  CHART  .  97 


iii 


LIST  OF  FIGURES 


Figure  Title  Page 

1  Impulse  Control  Guidance  Head  .  2 

2  Impulse  Control  Geometry  .  3 

3  Laser  Trajectory  Corrections  .  6 

4  Angle  Between  Longitudinal  Axis  and  Target  .  8 

5  Mil  Correction,  Angle  of  Attack,  Distance  to  Damp  versus 

Total  Impulse  . .  17 

6  Correction  versus  Thrust  Position  "L"  18 

7  Angle  of  Attack  versus  Time .  19 

8  Target  Array .  22 

9  Target  Miss  Statistics .  24 

10  Idealized  Thruster  Correction  .  28 

11  Alpha  versus  Beta  Body  Fixed  Coordinates .  29 

12  Alpha  versus  Beta  Fixed  Plane  Coordinates  .  30 

13  Average  Angle  of  Attack  at  Lase  Time  versus  Pitch 

Frequency . 31 

14  Angle  of  Attack  versus  Time,  Two  Pitch  Frequencies  .  34 

15  Average  Lase  Angle  of  Attack,  Mean  Point  of  Impact, 

and  CEP  versus  Pitch  Frequency  .  35 

16  Magnitude  of  Correction .  37 

17  Thruster  Correction  Error  . .  40 

18  Magnitude  of  the  Trim  Arm  versus  Roll  Rate .  46 

19  Average  Angle  of  Attack  at  Lase  -  Guided  Trajectories  ....  48 

20  Average  Number  of  Expended  Thrusters  versus  Spin  Rate  ....  49 

21  Total  Miss  Distance  (R  )  versus  Missile  Roll  Rate .  50 

»  L 

22  Total  Miss  Less  Gravity  Sag  (R  )  versus  Missile  Roll  Rate  .  .  51 

23  Pitch  Frequency  versus  Average  Launch  Velocity 

Baseline  Configuration  .  53 

24  Total  Miss  versus  Central  Dead  Zone  Size .  55 

25  Total  Miss  Less  Gravity  Sag  versus  Target  Velocity .  57 

26  Missile  Angle  of  Attack .  61 

27  Maximum  Angle  of  Attack  versus  Perturbation  Time .  63 

28  Sequenced  Pictures  Explosive  Blast  Wave  .  64 

29  Predicted  and  Measured  Angle  of  Attack .  66 


IV 


T'"  . . . 


LIST  OF  TABLES 


Table  Title  Page 

1  Baseline  Configuration  Physical  Characteristics  .  10 

2  Primary  Aerodynamic  Characteristics  .  10 

3  Performance  .  11 

4  Guidance  Summary  Statistics  .  26 

5  Accuracy  Data  versus  Pitch  Frequency . 32 

6  Magnitude  of  Correction  -  Nominal  Trajectory . 42 

7  Maximum  Expected  Correction  .  43 

8  Accuracy  Analysis  of  Two  Configurations  with 

Different  Expected  Correction  .  44 

9  30-Degree  Dive  Baseline  Configuration  for  Two 

Aircraft  Launch  Speeds . 54 

10  Baseline  Configuration  Accuracy  on  Standard  and 

Normal  Distribution  Arrays . 56 

11  Accuracy  for  Baseline  Configuration  for  Assumed 

Magnitude  of  Error  Sources . -58 


j 

i 

i 

t 


V 

(The  reverse  of  this  page  is  blank) 


SECTION  I 


INTRODUCTION 


During  the  past  2-1/2  years  a  study  has  been  made  to  determine  the 
feasibility  of  using  an  explosive  charge  to  give  terminal  correction  to  a 
fin  stabilized  rocket  using  a  fixed-body  seeker.  Figure  1  is  a  sketch  of 
the  type  guidance  head  being  considered.  Figure  2  is  a  sketch  of  the  rocket 
and  associated  four  quadrant  detecting  correction  method  to  be  employed. 

A  pulsed  laser  image  reflected  from  the  target  is  received  in  one  of  the 
four  quadrants  indicating  that  the  longitudinal  axis  of  the  missile  is  not 
in  line  with  the  target.  The  seeker  logic  then  fires  an  explosive  strip 
(action  time  approximately  10_<*  seconds)  in  the  image  quadrant.  The  ex¬ 
plosive  strips  are  expended  in  a  predetermined  sequence.  In  this  manner 
the  velocity  vector  (to  the  extent  that  the  angle  of  attack  is  small  when 
the  correction  is  made)  is  turned  in  line  with  the  target.  Part  of  the 
correction  comes  from  motion  of  the  eg  due  to  the  linear  momentum  explosive 
impulse,  and  the  remainder  comes  from  the  lift  force  generated  by  the  angle 
of  attack.  The  angle  of  attack  is  induced  primarily  by  the  moment  of  the 
impulse  which  is  forward  of  the  center  of  gravity. 

In  studying  the  guidance  system  it  is  necessary  to  simulate  the  motion 
of  the  rocket,  the  logic  and  correction  of  the  guidance  system,  and  all 
associated  errors.  Two  simulation  computer  programs  are  used  in  this  study. 
The  first  is  a  simple  closed  form  solution  to  the  pitching  motion  which  is 
used  in  analyzing  the  pitching  frequency,  angle  of  attack,  and  magnitude 
of  the  impulse  correction  for  the  high  force,  short  time  duration  impulses 
employed.  The  second  computer  program  is  a  complete  six-degree-of-freedom 
program  which  models  the  entire  guidance  system.  This  includes  physical  and 
aerodynamic  properties  of  the  rocket,  the  guidance  logic,  the  force  and 
moment  of  the  impulse  corrections,  atmospheric  conditions  and  target  motion. 

The  advantages  of  the  type  of  guidance  being  considered  are  that  it  is 
simple  and  inexpensive.  The  guidance  package  contains  no  gyros  or  other 
moving  parts. 

The  following  phenomena  can  adversely  affect  the  accuracy  of  this  body- 
fixed  system,  and  an  understanding  of  them  is  essential  in  developing  the 
system: 


(1)  Limited  correction  from  explosive  impulse  strips;  the  total 
correction  to  be  obtained  from  this  type  of  guidance  is  determined  by  the 
total  impulse  per  strip,  the  number  of  strips,  the  physical  and  aerodynamic 
properties  of  the  rocket,  and  the  type  of  trajectory  being  considered. 

(2)  Angle  of  attack  when  guidance  correction  made;  the  angle  of 
attack  of  the  rocket,  induced  by  impulse  correction  or  trim  misalignment, 
is  erroneously  interpreted  by  the  seeker  as  an  error  in  the  direction  of 
the  velocity  vector. 
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Figure  2.  Impulse  Control  Geometry 


(3)  Gravity  sag;  the  inability  for  the  logic  to  account  for  the 
normal  gravity  turn  of  the  rocket. 

(4)  Target  motion  and  wind;  the  system  does  not  take  into  account 
the  relative  motion  of  the  target  (no  lead  computing  capability),  and  this 
results  in  impacts  trailing  the  target. 

(5)  Disturbed  flow  due  to  the  blast  wave  from  the  impulse  explo¬ 
sive  strip;  if  the  disturbed  flow  from  the  explosive  strip  changes  the 
aerodynamic  forces  on  the  missile  for  a  long  enough  period  of  time,  an  un¬ 
predictable  correction  would  result. 

All  of  the  above  phenomena  will  be  considered  in  detail  in  this  report 
and  an  assessment  of  the  accuracy  of  the  system  will  be  given. 
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SECTION  II 


GENERAL  GUIDANCE  CONSIDERATIONS  AND  SIMULATION  MODELS 
1.  BASIC  GUIDANCE  SYSTEM 

The  guidance  corrections  are  initiated  by  the  seeker  determining  the 
off  axis  quadrant  of  a  reflected  laser  beam  from  the  designated  target. 

The  line  of  sight  laser  reflection  is  focused  by  the  lens  (Figure  1)  onto 
the  detector.  The  seeker  logic  determines  the  quadrant  of  the  focused 
beam  (Figure  2)  to  the  limits  of  the  dead  zone.  The  dead  zone  may  be  a 
small  area  at  the  center  of  the  field  of  view,  or  a  region  on  the  perpen¬ 
dicular  lines  separating  the  quadrants.  For  the  system  being  considered 
the  size  of  the  dead  zone  is  on  the  order  of  0.1  degree. 

The  seeker  receives  a  signal  each  time  the  laser  pulses.  The  primary 
laser  pulse  frequency  considered  in  this  report  is  20  per  second.  Upon 
receiving  a  signal  in  one  of  the  four  quadrants,  the  seeker  logic  will  desig¬ 
nate  an  explosive  strip  (thruster)  in  the  same  quadrant  to  fire  (Figure  2). 

In  determining  the  target  quadrant  the  seeker  has  determined  the  line  of 
sight  to  the  target  with  respect  to  the  longitudinal  axis  to  within  90 
degrees.  There  is  essentially  no  time  delay  between  reception  of  laser 
image  and  completion  of  thrust  action  (less  than  10-4  seconds),  and  this 
is  one  of  the  advantages  of  the  system.  The  logic  of  the  sequencing  of  the 
thrusters  may  be  varied.  Also  the  number,  the  total  impulse  and  the  location 
of  the  thrusters  may  be  varied.  These  variations  will  have  an  effect  on 
the  accuracy  and  will  be  discussed  in  detail  in  later  sections  of  this  report. 
The  amount  of  correction  which  can  be  obtained  from  a  configuration  of 
thrusters  is  basic  to  the  guidance,  and  is  discussed  in  Section  IV. 

The  guidance  being  considered  here  is  fixed  body.  This  means  there  are 
no  gyroscopes  or  any  other  moving  parts  to  identify  an  inertial  reference 
or  rocket  relative  motion.  The  guidance  system  can  only  determine  the 
angle  between  the  target  and  the  longitudinal  axis  (no  floating  head  as  in 
the  laser  guided  bomb).  It  is  desired  to  correct  the  direction  of  the 
velocity  vector  to  the  target,  but  by  the  nature  of  the  system  this  is  done 
only  to  the  extent  that  the  angle  of  attack  is  small  when  correction  is 
initiated.  Figure  3  illustrates  this  point  where  at  3.3  seconds  of  flight 
a  target  offset  50  feet  amounts  to  a  0.9  degree  correction.  If  the  angle 
of  attack  of  the  missile  is  any  sizeable  fraction  of  this,  then  unacceptable 
errors  result  (the  effect  on  accuracy  of  this  error  will  be  considered  later 
in  this  report) . 

Angle  of  attack  errors  can  be  induced  by  the  impulse  correction  itself. 

In  general  the  explosive  strip  will  be  located  forward  of  the  rocket  center 
of  gravity  as  shown  in  Figure  2.  As  a  result  of  the  explosive  impulse  the 
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Mil  Correction 


Figure  3.  Laser  Trajectory  Corrections 


missile  will  incur  an  angle  of  attack  which  will  oscillate  and  damp  ac¬ 
cording  to  the  physical  and  aerodynamic  properties  of  the  missile.  On  the 
following  laser  pulse  a  residual  angle  of  attack  will  be  interpreted  as  a 
target  direction  error  to  correct.  Since  we  are  attempting  to  correct  the 
velocity  vector  to  the  target,  accuracy  will  be  reduced.  The  relationship 
between  the  laser  pulse  frequency,  missile  pitch  frequency,  angle  of  attack, 
and  accuracy  will  be  discussed  in  Section  III. 

The  explosive  strip  when  detonated  will  also  cause  a  blast  wave  and 
associated  disturbed  flow  over  the  missile.  If  this  disturbance  affected 
the  missile  long  enough,  it  could  result  in  an  undesired  angle  of  attack 
which  could  reduce  accuracy.  A  test  to  determine  if  this  was  a  problem  was 
run  and  is  discussed  in  Section  VII. 

Another  source  of  angle  of  attack  error  is  trim  misalignment.  This 
results  in  an  angle  of  attack  in  the  trim  plane  rotating  with  the  spin 
frequency  of  the  missile.  The  relationship  between  the  static  trim  mis¬ 
alignment,  the  spin,  pitch  frequency,  and  the  resulting  trim  angle  will  be 
given  in  Section  V.  A  method  for  controlling  the  magnitude  of  the  trim  arm 
with  the  spin  and  thereby  reducing  the  guidance  error  is  also  given. 

Gravity  sag  is  a  source  of  error  in  this  system.  It  results  from  the 
guidance  system  correcting  to  the  line  of  sight  direction  to  the  target 
without  regard  to  the  normal  gravity  turn  of  the  missile.  This  is  depicted 
in  Figure  3,  where  at  guidance  initiation  (2.1  seconds)  the  velocity  vector 
is  pointed  126  feet  downrange  of  the  unguided  impact  point.  The  zero  angle- 
of-attack  missile  will  incorrectly  detect  a  target  at  the  unguided  impact 
point  as  being  about  0.5  degree  down.  Corrections  will  be  initiated  and 
will  result  in  the  missile  having  to  correct  back  to  the  target  in  the  latter 
part  of  the  trajectory.  The  missile  won't  be  able  to  correct  all  the  way 
back  to  the  target  and  will  impact  short.  Figure  4  shows  the  angle  between 
the  longitudinal  axis  and  the  target  as  a  function  of  space  slant  range  to 
the  target  for  a  nominal  rocket  configuration.  The  three  trajectories  are 
unguided  with  the  target  being  located  on,  50  feet  short  of,  and  50  feet 
beyond  the  unguided  impact.  Generally  the  gravity  sag  error  is  decreased 
by  increasing  the  total  impulse  per  correction.  This  will  be  shown  in  the 
six-degree-of-freedom  accuracy  study. 

The  guidance  system  being  considered  uses  pursuit  guidance  (i.e., 
correction  directly  toward  the  target).  This  results  in  the  missile  trailing 
the  target  if  there  is  relative  motion  between  the  missile  and  the  target. 

The  relative  motion  can  result  from  the  target  moving,  or  a  wind  causing  the 
rocket  to  move  with  respect  to  the  target,  or  both.  This  error  can  be  re¬ 
duced  by  increasing  the  total  impulse  per  thrust.  This  will  be  demonstrated 
in  Section  III. 


2.  NOMINAL  TRAJLCTORY,  BASELINE  ROCKET 


The  final  configurations  of  the  unguided  and  guided  version  of  the 
Advanced  Tactical  Rocket  (ATR)  have  not  been  determined  at  this  time.  The 
approximate  dimensions,  mass,  moments  of  inertia  and  aerodynamic  coefficients 
will  be  given  below.  These  data  are  based  on  preliminary  design  work  and 
aerodynamic  test  at  AEDC  (References  1  and  2) .  Most  of  the  analysis  in  this 
and  later  sections  will  deal  with  the  baseline  rocket  and  modifications  of 
this  configuration. 

The  rocket  being  considered  is  a  4- inch  diameter  rocket  which  burns  for 
about  2  seconds  giving  an  axial  impulse  of  approximately  5560  lb-sec.  The 
physical  characteristics  of  the  baseline  rocket  are  given  in  Table  1,  the 
primary  aerodynamic  properties  in  Table  2,  and  the  performance  for  two  dive 
angles  in  Table  3.  All  symbols  are  defined  following  the  tables.  For  the 
30-degree  dive,  450  knot,  6000-foot  aircraft  launch  conditions,  the  base¬ 
line  rocket  at  burnout  has  a  velocity  of  3295  ft/sec.  The  burnout  mass 
is  1.759  slugs,  the  eg  is  4.697  feet  from  the  base  and  the  transverse  mo¬ 
ment  of  inertia  is  4.997  slug-ft-^.  At  guidance  initiation  (2.1  seconds), 
the  velocity  is  3278  ft/sec,  the  altitude  3769  feet,  the  velocity  vector 
has  turned  over  from  the  initial  30  degrees  to  31.9  degrees  and  the  slant 
range  to  unguided  impact  has  shortened  from  approximately  11,300  feet 
initially  to  7030  feet. 
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TABLE  1.  BASELINE  CONFIGURATION  PHYSICAL  CHARACTERISTICS 


t 

(sec) 

M 

(slugs) 

I 

X 

(slug-ft2) 

I 

(slug-ft2) 

X 

eg 

(ft  base) 

-L 

(ft  eg) 

T 

(lbs) 

0 

2.0 

2.477 

1.759 

0.0643 

0.05425 

8.01 

4.997 

3.909 

4.6967 

2780 

2780  terminate 

2.1 

GUIDANCE  COMMENCE 

1.6 

TABLE  2.  PRIMARY  AERODYNAMIC  CHARACTERISTICS 


4ach 

C 

X 

C 

na 

C 

P 

(ft)  Base. 

C 

m 

a 

C 

mq 

C£6 

% 

2.0 

0.49 

11.00 

1.792 

-95.84 

-4108 

0.075 

-6.4 

2.5 

0.41 

9.396 

2.008 

-75.79 

-3540 

0.052 

-6.25 

2.75 

0.377 

8.995 

2.154 

-68.61 

-3427 

3.0 

0.344 

8.594 

2.30 

-61.79 

-3315 

0.035 

-6.0 

3.5 

0.31 

7.964 

2.617 

-49.69 

-3156 

0.015 

-3.9 

f  =  1 

c  See  Appendix  C  for  force  and  moment  equations 


TABLE  3 .  PERFORMANCE 


LIST  OF  SYMBOLS 


t  =  time  from  launch 

M  =  mass  of  rocket 

I  =  axial  moment  of  inertia 

x 

I  =  transverse  moment  of  inertia 

y 

X  =  distance  to  eg  from  base 

eg 

-L  =  forward  distance  to  the  thruster  from  the  eg 

T  =  axial  thrust 

Cx  =  axial  drag  (thrust  off) 

C.,  =  normal  force  coefficient  derivative  due  to  ot 

Na 

C  =  restoring  moment  coefficient  derivative  due  to  a 

ma  6 

C  =  damping  moment  coefficient  derivative  due  to  q 

Alt  =  '  altitude 


X 

SL 

V 

X 

Z 


down  range  distance 
slant  range 

velocity  of  the  rocket 

X  component  of  velocity  of  eg 

Z  component  of  velocity  of  eg 

spin  driving  moment  coefficient 

spin  damping  moment  coefficient 

fin  cant  parameter  (1.0  for  this  configuration) 


W1 


3.  LINEAR  PITCHING  AND  SWERVE  MOTION 

The  linear  equations  of  motion  for  a  rolling  missile  are  described  in 
many  publications  including  Reference  3.  Appendix  B  gives  a  simplified 
version  of  the  equations  for  angular  motion  for  a  rolling  fin  stabilized 
missile. a  For  the  slow  roll  rates  considered  (10  to  20  hertz),  the  nutation 
is  only  about  1  percent  greater  than  the  precession  frequency.  The  frequency 
for  pure  pitching  motion  will  be  between  these  two  frequencies.  Then  for 
the  more  complicated  case  of  motion  induced  by  side  impulse,  the  planar 
motion  solution  will  give  most  of  the  information  required  to  understand 
the  motion. 


For  the  impulse  guidance  being  considered  it  is  important  to  understand 
the  relationship  between  the  motion  of  the  rocket  after  side  impulse,  and 
all  of  the  parameters  affecting  the  motion.  The  equations  determining  the 
angular  and  swerving  motion  for  the  planar  motion  case  are  derived  in  Ap¬ 
pendix  A  (see  Figure  A-l).  In  the  derivation  of  these  equations  the  usual 
small  angle  and  magnitude  simplifications  are  made.  Steady  state  conditions 
(velocity,  spin,  etc)  are  also  assumed. 


The  quantities  of  importance  are  the  angle  of  attack  a.  the  pitch 

l  ^  ^  ) 

frequency  (oj  )  and  the  side  deflection  or  swerve  (<$ _) .  These  quantities 

1  5 

are  listed  in  the  following  equations: 


a 


(K 

K.ot  \ 

n  r»nc  t  + 

_L  _ 

a  t 

sin  0)^2 

U  bU  j  uj  -  L  a  ' 

T  12 

l  “l 

2u>i  y 

_ 

.. 

(A-23) 


a 

T 


At 


-cos  CJ^T 


sin  u^t 


+  1 


A  = 


(A-19) 


Tfte  equation  in  Appendix  B  omits  magnus  terms  and  terms  which  are  con¬ 
sidered  small  for  the  high  static  margin  missile  considered  in  Section  2. 

^In  computing  the  swerve  for  a  rolling  missile  an  average  thrust  is  com¬ 
puted.  Equation  (A-32)  Appendix  A.  The  linear  spin  equation  is  also  given 
in  Appendix  A. 
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CA- 20) 

(A- 15) 

(A- 16) a 

(A-10) 

CA— 12) 


.  -fv  _  r  xt  Kd 

=  ~~r~  cm  a  e  Tv —  sm  “i 

s  d  Na  p  4K  u,  1 

r  a  1 


T  +  2t  + 


(A-29) 


6^  =  Change  in  velocity  vector  direction  in  radians 

^  =  Missile  angle  of  attack  at  time  t^ 


=  Thrust  termination  time 


aIt  is  of  particular  interest  that  the  second  term  in  involves  C  alone 

P  mq 

without  Cm^.  This  suggests  that  it  may  be  possible  to  extract  Cmq  separately 
from  free  flight  range  test  to  the  extent  that  Cm  and  Fs  are  known. 
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■  Ti"e  after  Slde  thrust  termination,  t.,  =  0  when  t  ,  T 


The  missile  angle  of  attack  at  side  thrust  termination 


Missile  angle  of  attack  rate  at  side  thrust  termination  time 


=  Missile  mass 


=  Missile  velocity  (constant  in  above  equations) 
Air  density  (average) 


=  Missile  reference  area 


=  Missile  reference  diameter 


=  Dynamic  pressure  =  pV^/2 


Missile  axial  moment  of  inertia 


Missile  transverse  moment  of  inertia 


-  Average  constant  side  force  (See  Equation  (A-32) 


=  psd/2M 


=  Pitch  rate  (a  for  this  geometry) 


=  Angle  of  attack  rate 


=  25uS£fe{"  COnStant  in  the  ^Tferential  equation 


'  N0rmal  force  “efficient  derivative  due  to  a 


‘  DM,PinS  ’°ment  coefficient  derivative  due  to  q  (,  =  e) 


Cmd  =  Damping  moment  coefficient  derivative  due  to  i 


’  U  cgV°  aPPUCati0"  °f  Side  (negative 


the  pUcTf^ueJcTmun  Xoul 
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Section  III).  For  the  baseline  configuration  Equation  (A-15)  gives  an  av¬ 
erage  pitch  frequency  during  guidance  of  10.3  hertz.  The  angle  of  attack 
at  lase  time  will  be  some  fraction  of  the  maximum  angle  of  attack. 

Figure  5  is  a  plot  of  data  from  the  above  equations  for  the  baseline 
configuration.  The  maximum  angle  of  attack  (am)  is  obtained  from  Equation 
(A- 23)  and  is  plotted  as  it  varies  with  total  impulse.  It  is  desired  to 
keep  the  angle  as  small  as  feasible.  The  mil  change  in  the  velocity  vector, 
(6S  Equation  A-29) ,  as  it  varies  with  the  strip  impulse  is  also  plotted. 

This  value  should  be  large  to  obtain  large  side  corrections  (see  Section  IV) . 

The  distance  the  rocket  must  move  before  the  maximum  angle  of  attack 
damps  to  0.25  degree  is  also  plotted  in  Figure  5.  The  0.25  degree  was  se¬ 
lected  since  it  is  necessary  to  have  the  angle  of  attack  this  small  to  obtain 
sufficient  accuracy  (see  Section  III).  The  plot  was  made  from  Equation 
(A-23) .  For  reasonable  impulses  (4  lb-sec  and  above)  the  distance  to  damp 
to  0.25  degree  would  be  greater  than  700  feet.  This  means  that  we  cannot 
rely  on  damping  to  have  a  small  angle  of  attack  when  corrections  are  made. 

Figure  6  gives  the  mil  correction  and  maximum  angle  of  attack  for  the 
baseline  configuration  plotted  versus  the  distance  of  the  thrust  center 
forward  of  the  eg  location.  The  plot  shows  that  the  maximum  angle  of  attack 
is  extremely  small~(.005  degree)  if  the  side  thrust  is  located  at  the 
center  of  gravity  (L  =  0) .  The  mil  correction  for  this  thrust  location  is 
0.888  mil.  This  is  about  59  percent  of  the  nominal  value  (1.493  mils), 
corresponding  to  the  center  of  thrust  located  1.453  feet  (4.479  calibers) 
forward  of  the  eg.  Because  of  the  smaller  angle  of  attack  the  guidance 
system  is  more  accurate  for  the  thrust  center  close  to  the  eg  (this  is 
demonstrated  in  Section  III).  The  angle  of  attack,  the  mil  correction  and 
simplicity  of  design  (modular  concept)  must  be  considered  in  the  design. 

Figure  7  is  a  plot  of  the  angle  of  attack  versus  time  taken  from 
Equation  (A-23).  Also  plotted  on  the  same  curve  is  the  angle  of  attack  out¬ 
put  from  the  six-degree-of-freedom  program.  The  plot  shows  the  agreement 
is  very  good  for  the  first  1/2  cycle  after  guidance  initiation.  The  small 
difference  in  peak  value  is  due  to  a  slightly  smaller  constant  value  of  C 

mq 

being  used  in  the  planar  equations  than  the  Mach  dependent  C  used  in 

the  6-D  program.  The  angles  will  be  quite  different  after  2.15  seconds 
since  the  rocket  will  be  getting  thrust  of  various  phase  in  the  6-D  program 
while  the  rocket  in  the  linear  model  is  performing  pure  damped  pitching 
motion  from  one  impulse.  The  six-degree-of-freedom  motion  is  driven  by 
the  thrust  application  at  the  laser  pulse  frequency. 
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(Degrees) 


4.  SIX- DEGREE -OF -FREEDOM  GUIDANCE  SIMULATION  AND  ACCURACY  ANALYSIS  MODEL 

a.  Six-Degree-of-Freedom  Guidance  Simulation 

In  order  to  study  the  complex  guidance  problem  being  considered  it  is 
necessary  to  use  a  six-degree-of-freedom  computer  guidance  model.  The  basic 
model  was  taken  from  Reference  4.  This  computer  model  is  set  up  in  terms  of 
modules  which  can  be  used  independently.  The  modules  used  directly  from 
Reference  4  were  the  engine  module  (rocket  axial  thrust),  the  steady  wind 
and  air  data  modules,  translational  and  rotational  rigid  body  dynamics 
modules,  the  transformation  module  (body  fixed  J  inertial  coordinates),  and 
the  program  termination  and  miss  distance  calculation  module.  Integration 
of  the  state  variables  is  performed  using  an  Adams-Moulton  (fourth-order) 
integration  algorithm  with  Runga-Kutta  start  and  fixed  integration  step. 

The  equations  for  the  explosive  side  thrusting  forces  and  associated  moments 
(which  are  unique  to  the  system  being  considered  in  this  report)  and  the 
aerodynamic  forces  and  moments  were  written  as  separate  modules  in  this 
study.  The  accuracy  analysis  methods  and  equations  were  also  done  in  this 
study. 

A  general  model  of  the  forces  and  moments  generated  by  the  forward  po¬ 
sition  explosive  thrusters  has  been  set  up  and  programmed.  The  equations 
and  logic  of  this  fixed  body  model  are  given  in  Appendix  C.  The  model  is 
not  complicated  by  the  necessity  of  actuating  control  valves,  control  sur¬ 
faces,  or  other  complicating  mechanisms.  Upon  the  reception  of  a  simulated 
lazer  target  image,  error  signals  are  generated  in  the  fixed  body  reference 
system  according  to  Equation  (1) : 


-1 

e  =  Tan  R  /R 

q  z  x 

CD 

-1 

e  =  Tan  R  /R 

r  y  x 

e  =  Pitch  error  signal  of  target 

er  =  Yaw  error  signal  of  target 

R^  =  X  position  of  target  in  body  axis  system 

Ry  =  Y  position  of  target  in  body  axis  system 

R  =  Z  position  of  target  in  body  axis  system 

z 

ces  and  moments  are  then  applied  to  the  rocket  as  given  in  the  flow 
'gic  and  equations  of  Appendix  C. 


! 
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The  program  allows  for  a  variable  number  of  thrusters,  thrust  magnitude 
and  position  of  each  individual  thruster,  field  of  view,  dead  zone  size, 
and  sequencing  of  thrusters.  The  computer  program  has  been  expanded  to  in¬ 
clude  various  thruster  sequencing  and  other  options  not  included  in  Ap¬ 
pendix  C. 


In  order  to  obtain  adequate  accuracy  during  the  short  burn  time  of  the 
explosive  (order  of  10'4  seconds),  it  is  necessary  to  go  to  integration 
times  of  about  1/10  the  explosive  action  time.  To  reduce  the  computer  run 
ning  time  and  still  retain  accuracy  it  is  necessary  to  change  integration 
step  size  at  the  beginning  and  end  of  each  explosive  thruster  burn. 


The  body  axis  aerodynamic  forces  and  moments  are  also  listed  in  Appen¬ 
dix  C.  The  forces  include  the  axial  drag,  and  the  normal  forces  (due  to 
angle  of  attack).  The  moments  include  the  spin  driving  and  damping  moments, 
the  pitching  and  yawing  moments  due  to  angle  of  attack,  the  damping  moments 
due  to  body  rates  and  angle  of  attack  rates,  and  the  magnus  moments.  The 
axial  drag  and  spin  moment  coefficients  are  input  to  the  program  as  functions 
of  the  Mach  number  and  total  angle  of  attack  (T) .  All  other  coefficients 
are  input  to  the  program  as  functions  of  Mach  number  and  either  a  or  6  (angle 
of  attack  components,  see  Figure  C-6) 


a  =  Tan  [w/u] 

6  =  Tan_1[v/U] 

o  =  Tan"1 ^  [v2  +  w2 J  /  [u2  +  v2  +  w2  ] 

u,  v,  w  =  components  of  missile  velocity  in  x,  y,  z  body 
axis  system. 


b.  Accuracy  Analysis  Model 

In  evaluating  the  accuracy  of  a  configuration  of  the  system,  a  target 
position  is  simulated  in  the  ground  plane.  This  position  is  within  a  few 
hundred  feet  of  the  non-perturbed  unguided  impact  point  of  the  rocket.  In 
order  to  get  a  statistical  evaluation  of  the  accuracy,  it  is  necessary  to 

simulate  a  number  of  targets.  Figure  8  shows  a  normal  distribution  of  im¬ 

pact  points  for  the  unguided  version  of  the  rocket  (xg  and  yg  are  the 

inertial  coordinates  in  ground  plane) .  The  cross  marks  give  the  impact 
points  for  a  37.5  degree  dive  at  a  10,000-foot  slant  range  for  a  normal 
distribution.  The  standard  deviation  in  xg  and  yg  are  equal  to  7  mils  per¬ 
pendicular  to  the  trajectory  (impacts  in  xg  and  yg  are  assumed  independent). 

The  pattern  is  stretched  out  in  the  x  direction  since  the  7  mil  error  goes 
as  the  cosecant  of  the  dive  angle  in  the  ground  plane. 
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The  procedure  employed  in  evaluating  the  accuracy  is  to  use  the  unguided 
impact  points  as  targets  for  the  guided  version.  The  guided  version  would 
impact  at  the  center  of  the  array  if  guidance  were  not  activated.  In  order 
to  reduce  computer  time  and  better  analyze  accuracy  with  regard  to  target 
position,  the  13  circles  (Figure  8)  were  used  as  the  target  array  for  most 
of  the  accuracy  analysis  done.  One  of  the  target  circles  is  located  at  the 
center  of  the  pattern.  Eight  targets  are  50  feet  from  the  center  separated 
by  45  degrees  starting  with  a  downrange  position.  The  last  four  targets  are 
located  100  feet  from  the  center,  90  degrees  apart  starting  45  degrees  be¬ 
tween  x  and  y  .  It  is  shown  in  Section  VII  that  for  most  cases  the  difference 
e  e 

in  accuracy  between  the  symmetrical  target  array  (circles)  and  the  random 
normal  array  (marks)  is  small. 

The  accuracy  of  a  configuration  is  stated  in  terms  of  a  mean  point  of 

impact,  and  a  circular  error  probable  (CEP)  radius.  Figure  9  shows  the 

mean  point  of  impact  (with  respect  to  the  targets)  for  the  array  of  13 

targets  (circles  Figure  8).a  The  data  is  for  baseline  configuration,  a  30- 

degree  dive  from  6000  feet.  Ax  and  Ay  are  the  x  and  y  miss  distances 

6  6  6  6 

for  the  13  trajectories.  The  standard  deviations  of  the  impact  point  are 
calculated  by  the  following  formula. b 


Ax 


ep 


Y (AX  .  -  AX  )2 
/  v  epi  ep' 


N  -  1 


where 


___ 

B 


AY  .  -  AY  V 
ex  e ) 


N  -  1 


AX  .  ,  AY  .  = 

epi  ei  . 

AX  ,  AY  = 
ep  e 


the  difference  between  component  impact  and 
target  position  for  target  i 

the  mean  value  of  AX  .  and  AY  . 

epi  ei 


°Axep  anc*  °Aye  are  converted  to  a  CEP  according  to  Reference  5.  (See  Fig¬ 
ure  C-5) .  The  computer  program  prints  out  the  mean  point  of  impact,  the 
standard  deviation  of  the  impacts,  the  CEP,  and  other  data  of  interest  for 
the  total  number  of  targets  in. the  target  array.  The  same  data  can  be  ob¬ 
tained  for  any  combination  of  the  targets  in  the  array. 


aThe  AX  .  listed  means  the  miss  distance  perpendicular  to  the  trajectory. 

6p  1 

^For  small  N,  N-l  is  replaced  by  N. 
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Table  4  is  the  computer  printout  for  the  target  array  given  in  Figure  8. 
For  this  particular  case  the  target  array  was  broken  down  into  five  groups 
(information  on  miss  distance  for  long  as  compared  to  short  targets  can  be 
obtained  in  this  manner).  The  fourth  column  is  for  all  13  targets.  The 
first  row  is  the  missile  average  angle  of  attack  when  a  target  image  is 
received  (at  lase  time).a  The  standard  deviation  of  the  various  means  are 
listed  in  the  row  following  the  statistic. 


The  perpendicular  miss,  AX  mean,  is  the  miss  distance  perpendicular  to 

the  trajectory  and  is  obtained  by  multiplying  the  average  miss  in  the  hori¬ 
zontal  plane  by  the  sin  of  the  dive  angle.  The  y  mean  is  AYe.  The  radial 


miss  distance  is  the  distance  to  the  mean  point  of  impact  as  listed  in 
Figure  9.  The  data  in  Table  4  is  for  a  launch  of  30  degrees  from  6000  feet 
altitude. 


The  statistics  of  this  section  are  not  meant  to  portray  rigorous  sta¬ 
tistical  definitions.  The  methodology  is  set  up  to  evaluate  error  sources 
and  get  an  approximate  value  for  the  system  miss  distances. 


aAll  angles  are  in  degrees  and  distances  in  feet. 
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TABLE  4.  GUIDANCE  SUMMARY  STATISTICS 


Target  Sub  Arrays 


1-9 


Mean  Lase  Time  Angle  of  Attack  a £  (deg)  0.077 
Standard  Deviation  (deg)  0.013 


Mean  Perpendiculaj  Miss  Distance 
(Down  Range)  AXgp  (ft) 

Standard  Deviation  (ft) 


Standard  Deviation  (ft) 


Mean  Radial  Miss  Distance  R  (ft) 


Circular  Error  Probability  acep  (ft) 


Total  Miss  Rt  =  R  +  ocep  (ft) 


Mean  Dead  Zone  Occurrences 


Mean  Overexpenditures  of  Thrusters 


I 

6-9 

10-13 

0.081 

0.079 

0.008 

0.004 

-3.03 

-1.55 

1.24 

0.56 

1.28 

0.36 

0.66 

0.19 

3.29 

1.59 

1.10 

0.44 

4.39 

2.03 

19.0 

17.0 

0.25 

0 

kill 


1-13 


0.079 

0.010 


SECTION  III 


RELATION  BETWEEN  THE  LASER  PULSE  FREQUENCY  AND  PITCH  FREQUENCY 

As  mentioned  previously  in  this  report,  the  angle  of  attack  of  the  missile 
must  be  relatively  small  when  a  correction  is  made  since  the  angle  of  attack 
is  interpreted  as  an  error  in  the  velocity  vector  direction.  It  was  also 
shown  in  Section  II  that  there  was  not  sufficient  damping  between  laser  pulses 
to  reduce  the  angle  of  attack  to  acceptable  values.  The  method  used  to  over¬ 
come  this  problem  is  to  have  the  missile  pitch  frequency,  (essentially  the 
same  as  nutation  and  precession  frequencies  for  slow  roll  see  Section  II), 
slightly  larger  than  1/2  the  laser  pulse  frequency  (LPF).  Under  these  con¬ 
ditions  the  angle  of  attack  tends  to  be  small  when  a  correction  is  made,  and 
the  angle  of  attack  error  (due  to  pitch  frequency  being  above  1/2  the  LPF) 
is  opposite  the  correction  side  of  the  previous  correction.  The  latter 
tends  to  cause  out  of  phase  corrections  and  prevents  resonant  angle  of  attack 
buildup. 

The  above  can  be  qualitatively  understood  from  Figure  10. a  A  correction 
is  made  toward  the  target  where  the  initial  velocity  vector  is  V^.  The  angle 

of  attack  moves  to  and  swings  back  through  the  zero  angle  of  attack 

position,  (since  the  pitch  frequency  is  a  little  over  1/2  the  LPF).  The 
second  correction  is  made  when  the  angle  of  attack  is  «2£  and  the  correction 
F2  is  generally  in  the  opposite  direction  of  the  angular  velocity.  This 

tends  to  cause  the  random  error  corrections  to  be  out  of  phase  with  the  an¬ 
gular  velocity  and  prevents  a  resonant  buildup  of  angle  of  attack.  If  the 
pitch  frequency  were  less  than  1/2  LPF  a  resonant  buildup  of  the  angle  of 
attack  would  occur,  and  the  accuracy  would  be  greatly  reduced.  This  is 
clearly  demonstrated  from  the  6-D  output  given  below. 


Figure  13  is  a  plot  of  the  average  angle  of  attack  at  lase  (correction) 
time  for  15  variations  of  the  natural  pitch  frequency  of  baseline  configura¬ 
tion.  The  different  pitch  frequencies  were  obtained  by  varying  the  eg  loca¬ 
tion  of  this  configuration  (see  Table  5).  Each  of  the  12  points  in  Figure  13 
represents  the  results  of  running  trajectories  for  the  target  array  given  in 


aFigure  10  shows  planar  motion;  however,  for  the  high  impulse  system  being 
considered  the  six-degree-of-freedom  motion  is  nearly  planar.  6-D  plots  of 
the  angular  motion  (a  vs  8)  are  shown  in  Figures  11  and  12  (same  trajectory). 
Figure  11  is  in  the  body  fixed  coordinate  system  and  Figure  12  is  in  the 
fixed  plane.  Figure  12  shows  that  for  the  slow  roll,  and  out  of  phase  fast 
impulse  being  considered,  the  motion  is  nearly  planar  in  the  plane  of  the 
impulse. 
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Figure  10.  Idealized  Thruster  Correction 

angular  velocity 

maximum  angle  of  attack  on  first  swing 
velocity  prior  to  first  correction 

velocity  after  first  correction 

Angle  of  attack  when  second  correction  is 
at  second  laser  pulse 

second  thrust  force 
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Figure  11.  Alpha  Versus  Beta  Body  Fixed  Coordinates 
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Figure  12.  Alpha  Versus  Beta  Fixed  Plane  Coordinates 


TABLE  5.  ACCURACY  DATA  VERSUS  PITCH  FREQUENCY 


CG  Base 
(ft) 


.196 

.362 

4.445 

4.487 

4.5287 
4.57 
4.6967 
4.976 
5.0287 
5.2787 

5.5287 


Planar  Linear 
Theory  Data 

Six-Degrees-of-Freedom  Output 

“l 

(hertz) 

a 

m 

(deg) 

al 

(deg) 

R 

(ft) 

Hi 

8.83 

1.24 

1.71 

12.0 

19.8 

9.23 

1.63 

19.9 

16.5 

9.60 

1.16 

1.16 

20.8 

17.1 

9.78 

0.74 

18.7 

10.0 

9.86 

0.50 

15.5 

8.0 

9.95 

1.12 

0.24 

7.3 

3.0 

10.04 

0.042 

2.3 

0.5 

10.3 

0.079 

2.6' 

0.96 

10.85 

1.04 

0.22 

6.8 

2.5 

10.95 

0.21 

5.7 

1.8 

11.42 

0.30 

8.1 

2.7 

11.87 

0.37 

10.4 

4.4 

Figure  8  (13  targets).  The  plot  shows  that  above  the  natural  pitch  fre¬ 
quency  the  average  angle  of  attack  at  lase  is  related  to  the  angle  attained 
from  a  natural  swing  from  the  last  impulse. 


a,.  -  a  sin  aii»  =  a 
£2  m  cm 


sin 


2irfi 


(2) 


This  equation  gives  a  fairly  close  approximation  of  the  curve  (Figure  13) 
above  a  frequency  of  10  hertz.  Below  a  frequency  of  10  hertz  the  angle  of 
attack  builds  up  drastically.  This  is  due  to  the  bias  angle  of  attack  error 
being  in  the  direction  to  generally  cause  subsequent  impulses  to  be  in  phase 
with  the  angular  velocity.  This  creates  a  large  buildup  in  angle  of  attack. 
This  is  also  shown  in  Figure  14  where  the  upper  curve  is  for  the  frequency 
of  10.3  hertz  (above  resonance)  and  the  lower  curve  is  for  a  frequency  of 
9.78  hertz  in  the  resonant  region.  It  is  obvious  that  for  the  9.78  pitch 
frequency  the  angle  of  attack  is  building  up  and  is  up  to  about  2.0  degrees 
at  2.86  seconds  (the  plot  covers  only  about  1/3  the  guided  flight).  The 
plot  shows  the  average  angle  of  attack  at  lase  time  is  also  building  up. 

From  Table  5  it  is  seen  that  the  average  angle  of  attack  at  lase  for  the 
9.78  pitch  frequency  is  0.74  degree  compared  to  0.079  for  the  10.3  pitch 
frequency. 


The  effect  of  the  angle  of  attack  and  resonant  condition  on  guided  ac¬ 
curacy  are  shown  in  Table  5  and  Figure  15.  Table  5  lists  the  frequencies 
(o^)  and  maximum  angle  of  attack  (am)  from  the  planar  equations  (Section  II), 

and  some  of  the  statistics  from  the  6-D  output  summary  statistics  (mean 
angle  of  attack  at  lase  (a») ,  mean  point  of  impact  (R) ,  and  a  ). 

v  vCp 

The  distance  from  target  to  mean  point  of  impact  and  CEP  values  are  plotted 
versus  the  pitch  frequency  in  Figure  15.  The  angle  of  attack  is  again 
plotted  to  illustrate  the  relationship  between  resonance  angle  of  attack  and 
guided  accuracy  of  this  system. 
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Figure  14.  Angle  of  Attack  Versus  Time,  Two  Pitch  Frequencies 
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SECTION  IV 


APPROXIMATE  MAXIMUM  EXPECTED  CORRECTION 

1.  CORRECTION  ASSUMING  IMPULSE  EXACTLY  TOWARD  TARGET 

In  the  system  being  considered,  the  correction  toward  the  target  is  a 
result  of  the  thrust  impulses  from  the  explosive  strips  in  the  guidance 
head.  The  explosive  strips  are  limited  in  number  and  total  impulse  and  the 
maximum  expected  correction  toward  the  target  will  be  limited. 

The  angular  change  in  the  direction  of  the  velocity  vector  on  each  thrust 
will  result  from  the  rigid  body  impulse  and  the  lift  on  the  missile  due  to 
the  impulse  induced  angle  of  attack.  The  correction  will  be  a  function  of 
the  magnitude  of  the  impulse,  the  distance  from  the  rocket  eg  to  the  center 
of  the  impulse,  and  the  physical  and  aerodynamic  properites  of  the  rocket. 

The  equation  for  the  velocity  vector  angular  correction  is  derived  in  Sec¬ 
tion  II  and  Appendix  A  (Equation  A-29).  All  the  parameters  involved  in  ob¬ 
taining  the  angular  correction  (6S)  are  given  in  the  equation.  This  equation 
was  derived  using  a  planar  linear  theory  model,  however,  the  results  obtained 
from  using  the  equation  compare  favorably  with  the  data  from  the  six-degree- 
of- freedom  program. 

Equation  A-29  gives  the  angular  correction  to  the  velocity  vector,  but 
it  is  desired  to  get  the  distance  correction  on  the  ground.  In  order  to  get 
the  distance -correction  on  the  ground  it  is  necessary  to  consider  the  dis¬ 
tance  from  the  impact  when  the  correction  is  made.  Figure  16  is  a  pictorial 
view  of  the  missile  being  corrected  at  each  laser  impulse  as  it  approaches 
the  target.  The  approximate  distance  to  impact  at  the  i'th  lase  is 

Ri~(Ro  -  lV)  (3) 

where 

Rj  =  Approximate  distance  to  impact  at  the  i'th  lase 

i  =  i'th  lase.  i  =  0,  1,  2,  ...N 

Rq  =  Distance  to  impact  at  guidance  initiation 

V  =  Average  rocket  velocity  to  impact  determined  from  time  of 

flight  and  Rq 

T^  =  Period  between  laser  pulses 
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is  a  function  of  the  physica 
aerodynamic  properties  of  t 


The  ground  perpendicular  correction  due  to  the  i'th  lase  is 


=  <5  R. 


s  1 


(4) 


6  =  Angular  correction  to  velocity  vector  in  radians  (See 

s  Equation  A-29) . 


The  total  correction  is  the  sum  of  these  individual  corrections  through  the 
last  correction  before  impact. 


In  computing  the  total  maximum  correction  it  will  be  assumed  that  the 
guidance  head  has  a  central  dead  zone,  and  that  the  target  is  removed  enough 
from  the  unguided  impact  point  so  that  the  target  image  does  not  fall  in  the 
dead  zone.  Under  these  conditions  the  total  number  of  impulses  will  be  the 
number  of  thrusters  (N)  when  the  target  is  at  sufficient  range  so  that  the 
guided  flight  time  (Tf)  is  greater  than  (N-l)T^.  The  maximum  value  of  i 

(See  Figure  3)  for  this  case  is  (N-l).  For  shorter  ranges,  T^  <  (N-l)T^, 

the  maximum  value  of  i  is  the  integer  part  of  (T^/T^) .  The  maximum  value 

of  i  is  given  in  the  following  equations: 


=  Integer 
Part  of 


for 


~  <  (N-l) 

ll 


(5) 


M  =  (N-l), 


for 


>  (N-l) 

l 


where : 


M  =  Maximum  value  of  i  (See  Figure  3) 
T^  =  Guided  flight  time. 


The  average  velocity  is  given  in  terms  of  the  time  of  flight  by  the 
following  equation: 


V  =  R  /T- 
o  f 


(6) 
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The  total  correction  is  now  obtained  by  summing  Equation  (4)  from  zero 
to  M  where  M  is  given  by  Equation  (5)  and  V  by  Equation  (6).a 


2>i  -  CM  *  1)RC 


1  ML 

if 2  r- 


2.  EFFECT  OF  FOUR  QUADRANT  DETECTOR 

Since  the  target  is  detected  to  within  90  degrees,  the  correction  toward 
the  target  will  in  general  be  in  error.  This  is  demonstrated  in  Figure  17 
where  a  hypothetical  case  for  one  quadrant  is  shown.  The  correction  for  the 
target  image  shown  is  applied  at  the  angle  ef  rather  than  the  correct  angle 
e  .  Where  is  the  quadrant  angle  defining  the  direction  of  the  applied 

thrust  and  e  is  the  angle  defining  the  direction  to  the  target. 

For  better  values  for  correction  to  target,  the  error  angle  (e  -  Ef)  , 

tg  r 

of  the  applied  thrust  direction  should  be  taken  into  account.  The  thrust 
will  then  have  the  following  component  in  the  direction  toward  the  target. 

Ftg  =  F  C0S  (£tg  "  £f3  (8) 

Ftg  =  Component  of  thrust  toward  the  target 

-  e£  =  The  error  angle  of  the  applied  thrust 

The  thrust  will  also  have  a  component  perpendicular  to  the  target 
direction. 

F  =  F  sin  (e^  -  e-) 

P  1  tg  f  (9) 

Fp  =  Component  of  the  thrust  perpendicular  to  the  target  direction 

It  will  be  assumed  that  the  perpendicular  component  of  the  thrust  will 
average  to  zero  since  the  thrust  is  just  as  likely  to  be  on  one  side  of  the 
target  as  the  other.  It  is  necessary  to  find  the  expected  value  of  the 
error  angle  in  order  to  compute  the  expected  component  of  the  thrust  toward 
the  target. 


aIn  the  derivation  of  this  equation  the  fact  that  ^  \  =  M(M  +  l)j 2  was  used. 
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The  expected  value  of  the  error  angle  will  depend  on  the  geometry  and 
sequencings  of  the  thrusters  and  the  target  image  distribution.  It  is  rea¬ 
sonable  to  assume  that  the  angle  position  of  the  target  image  is  a  random 
variable  uniformly  distributed  over  the  quadrant.  In  order  to  simplify  the 
problem,  it  will  be  assumed  that  the  thrusters  also  have  a  random  uniform 
distribution  over  a  quadrant  (that  is,  no  account  will  be  taken  of  the  change 
in  the  thruster  distribution  with  thruster  depletion) . 

The  problem  and  solution  is  as  follows: 


FIND: 


"  er 

tg  f 


Expected  value  of  the  magnitude  of  the 
difference  between  the  thrust  and  target 
angle. 


Assume  a  uniform  distribution  of  e  and  e^.  Also  assume  e^g  is  independent 
of  e-,  then  the  probability  density  functions  for  e_  and  e  are 

t  ^  tg 

f  (e  )  =  1/90  ,  0  <  ef  <  90 


where 


f  (etg)  =  1/90  ,  0  <  ctg  <  90 

f  =  Probability  density  function. 


Assume  <(jetg  -  ef£> 


'  2\ 
S  -  £f>  )> 


Since  e^  is  independent  of  e^,  it  follows 


,(£tg  "  Ef 


o  -  <0  -2  O  <£f>  ♦  <4> 


Using  normal  procedure  in  evaluating  expected  value,  it  follows  that 


yu 

<4>  ■  <4>  ■/  S* 


dc  -m 

tg  3 


41 


and 


90 

<«>  ■  <■»  ■  / 


-I&de 

90  QGtg 
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(14) 


Substituting  in  Equation  (12)  gives 


(£tg  - 


*  37.74  degrees 


(15) 


then  the  expected  correction  toward  the  target  is 


6  cos  37.74  =  6  (.792) 

rp  rp  ' 


(16) 


3.  MAXIMUM  EXPECTED  CORRECTION  FOR  NOMINAL  CONFIGURATIONS 

The  utilization  of  the  maximum  correction  equations  will  be  demonstrated 
in  this  section.  All  the  rocket  systems  which  have  been  considered  give  the 
velocity  angle  correction  (Sg)  between  1  and  2  milliradius.  Then  from 

Equations  (7)  and  (16)  the  magnitude  of  the  perpendicular  correction  will  be 
computed  for  a  one  milliradian  system  (6g  =  0.001  radians).  Table  6  shows 

the  corrections  for  the  nominal  trajectories  listed  in  Table  3. 

TABLE  6.  MAGNITUDE  OF  CORRECTION  -  NOMINAL  TRAJECTORY 


Dive  Angle 

6 

6 

/6  \ 

s 

rp 

\  rp/ 

(degrees) 

(radians) 

(feet) 

(feet) 

30 

0.001 

159 

126 

37.5 

0.001 

96 

76 

An  example  of  using  the  maximum  correction  equations  is  given  in  the 
following  analysis.  The  planar  motion  equations  and  the  maximum  correction 
Equation  (16)  were  used  to  evaluate  two  guidance  configurations  for  the 
37.5-degree  dive  listed  in  Table  3.  The  configurations  were  the  same  except 
for  the  total  impulse  of  the  side  thrusters  (Ig)  and  the  center  of  thrust 
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location.  The  two  configurations  are  listed  in  Table  7.  The  table  shows 
that  the  velocity  vector  correction  (6g)  is  higher  on  Configuration  B  giving 

an  expected  maximum  correction  of  108  feet  compared  to  84  feet  on  Con¬ 
figuration  A.a 


TABLE  7.  MAXIMUM  EXPECTED  CORRECTION 


Planar 

Linear  Equations 

Expected 

Correction 

Equation 

I 

s 

L 

“l 

am 

6 

s 

<!r> 

Configuration 

(lb  sec) 

(feet) 

(hertz) 

(degrees) 

(radians) 

(feet) 

A 

4 

1.174 

10.86 

0.66 

0.00111 

84 

B 

4.4 

1.4 

10.86 

0.951 

0.00142 

108 

The  six-degree-of-freedom  guidance  model  was  used  to  evaluate  the  accuracy 
of  these  two  configurations.  The  target  array  used  was  that  given  in  Figure  8. 
This  array  has  five  targets  within  50  feet  and  four  targets  at  100  feet.  The 
two  sets  were  analyzed  separately  and  together.  These  results  are  shown  in 
Table  8.  Table  8  shows  that  for  configuration  A  the  CEP  values  get  quite 
large  (about  11  feet]  for  targets  at  100  feet.  This  is  readily  understood 
since  the  maximum  expected  correction  is  calculated  to  be  84  feet.  In  con¬ 
figuration  B  the  maximum  expected  correction  is  brought  up  to  108  feet  by 
increasing  the  total  impulse  and  moving  the  center  of  thrust  forward.  This 
enables  the  rocket  to  reach  the  100- foot  targets  and  the  CEP  comes  down  to 
a  low  value  (about  1.3  feet).  The  analysis  indicated  that  the  planar  equa¬ 
tions  and  the  simple  maximum  correction  equations  can  be  used  in  design  of 
the  impulse  guided  system. 


aThe  pitch  frequencies  (co^)  are  the  same  where  the  maximum  angle  of 
attack  (am)  is  higher  on  configuration  B. 


TABLE  8.  ACCURACY  ANALYSIS  OF  TWO  CONFIGURATIONS 
WITH  DIFFERENT  EXPECTED  CORRECTION 


Con¬ 

figuration 


<6  \ 

\  rp/ 

(ft) 

Target  Set 

AX 

ep 

(ft) 

AY 

e 

(ft) 

R 

(ft) 

oAX 

ep 

(ft) 

oAY 

e 

(ft) 

acep 

(ft) 

50-ft  targets 

-2.4 

0.2 

2.4 

0.5 

0.5 

0.6 

84 

100-ft  targets 

-1.7 

1.8 

2.5 

7.6 

11.4 

11.2 

All  targets 

-2.2 

0.7 

2.3 

4.4 

6.7 

6.5 

50-ft  targets 

-2.7 

-0.3 

2.7 

1.4 

0.5 

1.1 

108 

100-ft  targets 

-2.0 

0.3 

2.0 

0.5 

1.8 

1.3 

All  targets 

-2.5 

-0.1 

2.5 

1.2 

1.1 

1.4 

SECTION  V 


TRIM  MISALIGNMENT  EFFECT 

A  misalignment  of  the  body  components  or  tail  fins  will  cause  aerodynamic 
forces  and  moments  on  the  rocket  which  results  in  the  rocket  trimming  out  at 
some  angle  of  attack  greater  than  zero.  For  a  spinning  rocket  this  trim 
misalignment  causes  the  rockets  longitudinal  axis  to  perform  a  lunar  type 
motion  about  the  velocity  vector  at  an  angle  of  attack  (trim  angle) .  For 
the  guidance  system  being  considered  this  will  result  in  error  correction 
and  premature  expenditure  of  thrusters  giving  poor  accuracy.  This  phenomenon 
was  investigated  and  a  method  of  controlling  trim  misalignment  accuracy  error 
was  found. 

The  magnitude  of  the  dynamic  trim  arm  is  primarily  a  function  of  the 
trim  moment,  the  physical  and  aerodynamic  characteristics  of  the  missile,  and 
the  missile  spin.  The  magnitude  of  the  dynamic  trim  arm  is  given  in  Equation 
(B-4).  This  equation  was  derived  in  Appendix  B,  where  only  a  trim  moment 
is  considereda,  and  other  simplifying  assumptions  were  made. 


(B-4) 


6m  =  Static  trim  misalignment  (the  trim  angle  at  zero  spin  rate) 

p  =  Spin  rate 

^1,2  _  Precession  or  nutation  rate 

for  the  relatively  low  spin  rate  being  used 

*x  *  V  =  Axial  and  transverse  moments  of  inertia 


Figure  18  was  generated  from  Equation  (B-4)  and  is  a  plot  of  the  trim  arm 
versus  the  missile  spin  rate.*5  For  the  relatively  low  spin  rate  being  used. 


aNo  loss  in  generality  in  regard  to  the  trim  angle  of  attack  is  encountered 
if  the  trim  force  term  is  omitted. 

^All  data  generated  in  this  section  are  for  the  baseline  configuration  given 
in  Section  II.  The  spin  was  varied  by  changing  the  spin  cant  parameter,  f£. 
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Computed  from  Equation  B-4 

Average  Angle  of  Attack  from  Unguided  6-D  Program 


Figure  18.  Magnitude  of  the  Trim  Arm  versus  Roll  Rate 


$1  2  is  essentially  equal  to  the  one-dimensional  pitch  frequency  (uii).  The 
» 

curves  in  Figure  18  were  computed  for  the  30-degree  trajectory  given  in 
Table  3.  The  average  conditions  during  the  guidance  phase  were  used  in  com¬ 
puting  the  pitch  frequency  by  Equation  (A-15).  This  value  for  uj  (10.29  hertz) 
was  used  in  Equation  (B-4) .  The  same  trim  moment  used  in  generating  Fig¬ 
ure  18  was  used  in  the  6-D  program  on  the  unguided  trajectories.  The  average 
angles  of  attack  obtained  are  also  plotted  in  Figure  18  and  are  in  excellent 
agreement  with  Equation  (B-4).  The  curve  shows  the  steep  buildup  in  angle 
of  attack  as  the  spin  frequency  approaches  the  pitch  frequency  (10.29  hertz). 
Curves  are  drawn  for  three  static  trim  misalignments  (0.286,  0.127  and  .057 
degrees).  A  non-rolling  missile  trims  at  the  static  trim  misalignment  angle. 

Guided  trajectories  with  varying  spin  and  trim  misalignment  were  run 
using  the  target  array  given  in  Figure  8.  Figure  19  is  a  plot  of  the  average 
angle  of  attack  at  lase  time  during  the  guided  phase  for  these  trajectories. 

The  curves  show  that  the  angle  of  attack  varies  in  the  same  manner  as  that 
given  by  Equation  (B-4)  (See  Figure  18). 

Also  related  to  the  trim  induced  angle  of  attack  is  the  depletion  of 
thrusters  in  the  trim  angle  of  attack  plane.  Figure  20  shows  the  average 
number  of  times  a  thrust  was  attempted  from  a  quadrant  in  which  the  thrusters 
had  been  expended.  These  curves  are  drawn  from  the  same  guided  trajectories 
mentioned  above.  The  curves  follow  the  same  trend  as  the  angle  of  attack 
curves  going  to  small  values  at  higher  roll  rates.  It  would  be  expected 
that  the  configurations  having  high  average  angles  of  attack,  and  related 
over  expended  quadrant  thrusters,  would  have  high  miss  distances. 

Figure  21  shows  the  total  miss  distance  follows  the  same  trend  as  the 
angle  of  attack  and  thruster  depletion  in  regard  to  missile  roll  rate. 

Figure  21  consists  of  plots  of  the  total  miss  distance  (mean  point  of  impact 
plus  the  CEP)  plotted  versus  the  missile  roll  rate.  The  total  miss  is  also 
plotted  for  the  missile  with  no  trim  misalignment. 

Figure  22  is  the  same  as  Figure  20  with  the  2.5-foot  gravity  sag  of  the 
nonperturbed  missile3  (missile  with  no  trim  misalignment)  taken  out.  The 
total  miss  distance  for  the  low  spin  (12.25  hertz)  and  high  trim  misalignment 
(0.29  degree)  is  22.5  feet  but  the  plot  shows  that  the  total  miss  for  any 
trim  misalignment  used  is  less  than  2  feet  for  spins  above  19.5  hertz.  This 
analysis  indicated  that  miss  distances  due  to  trim  misalignment  can  be  re¬ 
duced  to  acceptable  values  by  increasing  the  roll  rate  to  above  1.9  times 
the  pitch  frequency. 


aIn  using  this  type  of  guidance  system  the  laser  would  be  aimed  high  to 
eliminate  part  of  the  normal  gravity  sag  error. 


Figure  20.  Average  Number  of  Expended  Thrusters  versus  Roll  Rate 
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SECTION  VI 


OTHER  GUIDANCE  CONSIDERATION  AND  ACCURACY  ERROR  SOURCES 
1.  PARAMETERS  AFFECTING  ACCURACY 

Several  things  that  affect  the  accuracy  of  this  type  of  system  will  be 
considered  in  this  section.  These  are  the  gravity  sag,  the  aircraft  launch 
velocity,  the  dead  zone  size,  target  motion  and  wind. 

The  gravity  sag  is  the  target  miss  due  to  the  natural  gravity  turn  angle 
of  the  missile.  The  velocity  vector  of  the  missile  at  guidance  initiation 
will  be  pointed  downrange  of  the  missile  impact  point.  The  logic  of  the 
pursuit  type  guidance  system  being  used  will  incorrectly  call  for  down  cor¬ 
rections  for  a  target  at  the  unguided  impact  point.  This  will  result  in  the 
missile  being  short  of  the  target  in  the  latter  phases  of  the  trajectory  and 
unable  to  make  up  the  miss.  There  are  a  number  of  parameters  that  will  affect 
the  magnitude  of  the  gravity  sag  error.  Some  of  these  are  the  velocity  of 
the  missile  at  burnout,  the  angle  and  slant  range  to  the  target,  and  the 
number  and  magnitude  of  the  thrusters.  The  gravity  sag  for  the  baseline 
configuration  for  a  30-degree  dive  from  6000  feet  is  2.5  feet. 

The  aircraft  launch  velocity  can  affect  the  accuracy  of  the  missile  by 
changing  the  pitch  frequency  and  by  changing  the  gravity  sag.  The  pitch 
frequency  is  given  by  Equation  (A-15).  If  terms  of  small  magnitude  are  ne¬ 
glected  the  pitch  frequency  reduces  to 


ujaV  tl  pC  (M) 
1  f  ma 


The  pitch  frequency  (u^)  is  seen  to  depend  on  the  velocity  linearly 

and  Cma<  Cma  decreases  with  velocity  (in  the  Mach  range  of  interest)  but 

V  is  dominant  in  Equation  (17).  From  Section  III  it  was  shown  that  the 
pitch  frequency  must  not  drop  below  1/2  the  laser  pulse  frequency  to  avoid 
angle  of  attack  buildup  and  loss  of  accuracy.  The  rocket  must  be  designed 
so  that  normal  variation  in  aircraft  launch  does  not  change  the  pitch  fre¬ 
quency  to  less  than  1/2  the  laser  pulse  frequency.  Figure  23  is  a  plot  of 
the  pitch  frequency  versus  the  aircraft  launch  velocity  (this  plot  was 
constructed  using  Equation  (17)  and  unguided  trajectory  data).  The  value 
of  the  frequency  is  nearly  constant  from  the  initiation  of  guidance  until 

impact.  This  occurs  since  p  and  C  are  increasing  during  descent  as  V 

ma 

decreased.  The  plot  shows  that  the  pitch  frequency  goes  above  10  hertz 
at  about  350  knots  for  the  baseline  configuration.  Table  9  shows  the  accu¬ 
racy  for  a  300-knot  aircraft  launch  (pitch  frequency  below  10  hertz),  and 
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the  450-knot  aircraft  launch  (pitch  frequency  above  10  hertz).  The  table 
shows  that  the  lower  aircraft  launch  results  in  resonant  angle  of  attack 
buildup  and  poor  accuracy.  This  particular  configuration  (baseline  config¬ 
uration)  is  designed  for  aircraft  launches  above  400  knots.  If  it  were 
desired  to  design  the  rocket  for  lower  aircraft  launches  the  pitch  frequency 
would  be  changed  by  a  eg  change,  a  fin  change,  or  other  changes. 


TABLE  9.  30-DEGREE  DIVE  BASELINE  CONFIGURATION 
FOR  TWO  AIRCRAFT  LAUNCH  SPEEDS 


Aircraft  Launch  Speed 

Mean  Angle  at  Lase  a ^ 

R 

CEP 

(Knots) 

(Degrees) 

(Feet) 

(Feet) 

300 

0.28 

16.5 

7.8 

450 

0.079 

2.6 

0.96 

The  dead  zone  is  a  region  on  the  image  plane  where  the  reflected  laser 
beam  does  not  result  in  a  signal  to  fire  a  thruster.  The  design  can  be 
varied  so  that  the  dead  zone  is  in  the  form  of  a  perpendicular  cross  marking 
the  boundaries  of  the  quadrants  or  in  a  central  region.3  An  analysis  was 
made  to  determine  the  effect  of  a  central  dead  zone  size  on  the  missile 
accuracy.  Figure  24  shows  the  results  of  this  analysis.  Figure  24  shows 
that  for  the  slant  range  and  dive  angle  considered  (11,300  feet  and  30  degrees) 
the  optimum  dead  zone  size  is  about  0.13  degree.  If  the  dead  zone  is  too 
large,  corrections  are  not  made  for  targets  sufficiently  far  from  the  impact 
point  until  the  missile  is  too  close  to  the  target.  For  the  dead  zone  too 
small,  an  excessive  error  is  picked  up  because  of  gravity  sag  effect  and, 
or  depletion  of  thrusters. 


Target  motion  and  wind  are  similar  in  causing  accuracy  errors  in  pursuit 
type  guidance  systems.  The  error  is  due  to  the  missile  continually  cor¬ 
recting  directly  toward  the  target.  An  additional  error  comes  about  if,  due 
to  target  motion  or  wind,  the  final  position  of  the  target  relative  to  the 
initial  position  is  beyond  the  maximum  expected  correction.  Several  moving 
target  arrays  were  simulated  on  the  computer.  The  standard  target  array  was 
used  (Figure  8),  with  the  target  motion  beginning  at  guidance  initiation 
(2.1  seconds).  The  target  motion  was  perpendicular^  to  the  trajectory  to 
the  right  and  varied  from  0  to  40  ft/sec.  This  would  correspond  to  an  air 


aThe  configuration  used  in  doing  most  of  the  error  analysis  in  this  report 
used  a  cross  dead  zone  of  0.05  degree. 

l- 

It  should  be  noted  that  a  wind  or  target  motion  perpendicular  to  the  tra¬ 
jectory  gives  about  50  percent  greater  error  than  downrange  motion. 
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Total  Miss  R  (feet) 
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mass  wind  with  the  same  speed  moving  in  the  opposite  direction.  The  total 
miss  distance  less  gravity  sag  (R  +  CEP  -  GS)  is  plotted  versus  the  target 
speed  in  Figure  25. a  The  plot  shows  that  this  error  is  large  and  is  up  to 
17  feet  for  a  30  ft/sec  target  speed.  The  plot  shows  that  the  higher  impulse 
reduces  the  error  some.  It  should  be  noted  that  it  is  common  operational 
practice  (laser  guided  bomb)  to  eliminate  part  of  this  error  by  leading  the 
target.  This  correcting  procedure  would  be  easier  to  accomplish  for  a  moving 
target  than  for  a  wind. 

2.  ESTIMATE  OF  THE  ACCURACY  OF  THE  SYSTEM 

Throughout  this  report  the  standard  target  array,  (the  solid  circles  in 
Figure  8),  has  been  used  when  computing  accuracy.  A  more  realistic  target 
array  would  be  that  given  by  the  distribution  representing  the  accuracy  of 
the  unguided  rocket.  Let  us  assume  that  X  and  Y  (unguided  rocket  impact 
points)  are  normally  distributed  random  variables  with  standard  deviations 
of  7  mils.  This  would  give  the  distribution  represented  by  the  cross  marks 
in  Figure  8  for  a  30-degree  dive  at  a  slant  range  of  11,000  feet.  For  the 
baseline  configuration  this  gives  about  8  percent  of  the  targets  beyond  the 
maximum  expected  correction.^  The  accuracy  of  the  baseline  configuration  was 
computed  for  stationary  targets  using  the  normal  distribution  but  not  in¬ 
cluding  the  8  percent  of  the  targets  beyond  the  maximum  expected  correction. 

A  comparison  of  the  accuracy  statistics  of  the  standard  and  the  modified 
normal  distributions  are  given  in  Table  10. 


TABLE  10.  BASELINE  CONFIGURATION  ACCURACY  ON  STANDARD  AND  NORMAL 

DISTRIBUTION  ARRAYS 


Array 

AX 

AY 

R 

CEP 

Standard 

-2.5 

0.86 

2.6 

0.96 

Normal 

Distribution 

(Modified) 

-2.5 

0.14 

2.5 

1.6 

Table  10  shows  very  little  difference  between  the  two  arrays  except  the  CEP 
is  slightly  higher  for  the  normally  distributed  array.  The  standard  target 
array  will  be  used  in  the  following  approximate  accuracy  analysis. 


aThe  nominal  gravity  sag  of  2.5  feet  was  subtracted  in  constructing  the 
curve . 

^This  may  be  computed  approximately  by  the  following  formula. 

Pj,  =  exp  2  where  P^  is  the  percent  of  target  beyond  the 

maximum  expected  correction,  a  =  =  standard  deviation  error  in  x,  and 

y,  and  =  maximum  expected  correction  (see  Section  IV) . 


Baseline  Configuration  30  Degree  Dive 
Target  Motion  Perpendicular  to  Trajectory 


Gravity  Sag  versus  Target  Velocity 


The  accuracy  of  the  system  was  very  roughly  estimated  here  using  some¬ 
what  arbitrary  assumption  on  the  magnitude  of  the  error  sources.  The  error 
sources  considered  are  the  gravity  sag,  trim  misalignment,  target  motion, 
and  wind.  The  errors  were  assumed  to  be  independent.  It  was  assumed  that 
the  laser  beam  would  be  aimed  in  an  attempt  to  reduce  the  error  due  to  gravity 
sag,  target  motion,  and  wind.  The  accuracy  study  was  run  for  the  baseline 
configuration  launched  from  6000  feet  altitude,  at  an  angle  of  30  degrees, 
from  an  aircraft  moving  at  450  knots. 

A  standard  deviation  for  the  gravity  sag  error  was  taken  as  50  percent 
of  the  nominal  gravity  sag  or  1.25  feet.  A  standard  deviation  error  in  trim 
misalignment  was  taken  as  a  static  trim  of  0.143  degree  for  a  roll  rate  of 
19  hertz.  From  Figure  22  this  gives  an  error  of  1.2  feet,  and  this  was  re¬ 
solved  into  a  and  a  by  multiplying  by  cos  45  degrees.3  A  standard 

deviation  error  for  target  motion  was  taken  as  a  target  moving  at  20  ft/sec 

which  was  50  percent  corrected.  The  20  ft/sec  was  resolved  into  a  cross 

range  and  downrange  component  of  14.14  ft/sec.  From  Figure  24  this  gives 

a  target  motion  error  of  8.2  feet  and  if  this  is  corrected  by  50  percent 

then  the  a  error  is  4.1  feet  for  lateral  target  motion.  The  error  in  the 
ym 

downrange  direction  is  about  one  half  this  for  a  30-degree  dive  or  axm  is 
2.05  feet. 

A  standard  deviation  error  in  wind  was  taken  at  30  percent  corrected 

10  mph  (14.67  ft/sec)  in  a  direction  45  degrees  off  the  downrange  direction. 

The  components  of  wind  in  the  cross  range  and  downrange  directions  are  10.37 

ft/sec.  From  Figure  24  this  gives  a  wind  error  of  6.6  feet.  Then 

a  =  0.7  (6.6)  *  4.62  feet  and  o  is  one  half  this  or  2.31  feet.  These 
yw  1  xw 

values  are  listed  in  Table  11.  Table  11  also  gives  the  root  mean  square 
value  of  ox  and  o^  and  a  CEP  error  is  computed  according  to  Reference  5. 

TABLE  11.  ACCURACY  FOR  BASELINE  CONFIGURATION  FOR  ASSUMED  MAGNITUDE 

OF  ERROR  SOURCES 


0 

a 

X 

y 

Gravity  sag  (50  percent  corrected) 

1.25 

Trim  misalignment  static  trim  0.143  degree, 
roll  rate  19  hertz 

0.85 

0.85 

Target  motion  20  ft/sec  50  percent  corrected 

2.05 

4.1 

Wind  10  mph  30  percent  corrected 

2.31 

4.62 

3.44  6.23 

i 

CEP  5.63 


a  is  taken  as  downrange  and  o  is  cross  range  standard  deviation, 
x  y 

subscript  t  stands  for  trim,  g  gravity  sag,  m  target  motion,  and  w  wind. 


The 
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The  error  analysis  given  in  Table  11  was  don.  for  theWrge 

ure  8.  This  array  roughly  yepresen  targets  beyond  the  maximum 

(about  rp=  of S error 

the  target  was  illuminated  without  jitter. 
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SECTION  VII 


DISTURBED  FLOW  EFFECT  AND  TEST 


For  the  guidance  system  being  considered  it  is  important  that  the  ex¬ 
plosion  from  the  thruster  does  not  adversely  affect  the  aerodynamic  charac¬ 
teristics  of  the  missile.  A  test  was  run  at  AEDC  (Reference  2)  to  determine 
if  the  disturbance  caused  by  the  explosive  strip  reacts  with  the  flow  to 
change  the  aerodynamics  of  the  missile  causing  changes  in  the  angle  of  attack 
or  pitch  frequency.  The  time  that  the  blast  wave  (or  disturbed  flow)  remains 
over  the  missile  is  important  in  determining  its  effect. 

The  solution  to  the  linear  equation  (Appendix  A)  was  used  to  estimate 
the  effect  of  disturbed  flow  over  the  missile  as  a  function  of  the  duration 
of  the  disturbance  (this  work  was  accomplished  by  the  author  and  reported  in 
a  memo  in  September  1975).  If  the  impulse  time  of  the  thruster  is  very 
short  and  the  relative  magnitude  of  physical  properties  and  aerodynamics  of 
the  missile  are  considered,  then  Equations  (A-19)  and  (A-20)  reduce  to 
Equation  (18) 


a  ~ 


FLT 

I 


(18) 


and  the  angle,  of  attack  and  angular  rate  of  the  missile  at  the  end  of  burning 
does  not  depend  on  the  aerodynamics  of  the  missile  at  all,  but  only  on  the 
thrust,  thrust  time,  moment  distance,  and  missile  transverse  moment  of  inertia. 
The  thrust  time  being  considered  for  the  explosive  impulse  ( — 1 0 ~ 4  seconds) 
is  short  enough  so  that  the  above  conditions  prevail.  After  the  explosive 
impulse  the  disturbance  due  to  the  explosive  shock  wave  and  associated  per¬ 
turbed  flow  will  affect  the  missile  aerodynamics  in  an  unknown  manner  for  an 
unknown  period  of  time.  Equations  (A-13)  through  (A-23)  were  used  in  an 
iterative  manner  to  determine  the  type  of  motion  obtained  for  a  given  pertur¬ 
bation  as  a  function  of  time.  It  was  found  that  changing  the  damping  co¬ 
efficient,  C  ,  had  little  effect  on  the  maximum  angle  of  attack  the  missile 
mq  6 

obtained.  The  maximum  angle  of  attack  was  a  strong  function  of  the  variation 

in  C  and  the  time  C  was  perturbed, 
mot  ma  r 


For  the  purpose  of  analysis  was  reduced  by  a  factor  of  (0.01)  and 

Equations  (A-13)  through  (A-23)  were  used  to  compute  the  value  of  ^  and 

for  various  perturbing  times.  Equation  (A-13)  was  used  to  continue  the 

computation  of  the  angle  of  attack  after  the  perturbation  was  terminated. 
Figure  26  shows  the  damped  oscillation  for  missile  characteristics  similar  to 
the  baseline  configuration.  The  maximum  angle  of  attack  obtained  is  about 
1.18  degrees.  The  other  curve  in  Figure  26  gives  the  missile  motion  if  Cma 

is  reduced  by  a  factor  of  0.01  for  a  period  of  20  milliseconds  after  thrust 
termination.  The  maximum  angle  of  attack  obtained  here  is  1.91  degrees. 
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-1.2 


Figure  26.  Missile  Angle  of  Attack 


The  simulations  show  that  the  effect  on  the  missile  motion  is  a  strong 
function  of  the  time  the  perturbing  influence  acts.  Figure  27  shows  the 
missile's  maximum  angle  of  attack  as  a  function  of  the  perturbation  time 
where  the  perturbation  is  a  reduction  in  Cma  by  a  factor  of  0.01.  This 

curve  shows  that  it  takes  a  perturbing  time  of  about  2  milliseconds  before 
any  change  in  maximum  angle  of  attack  is  noted.  After  this  time  the  effect 
multiplies  rapidly.  Although  here  the  perturbing  influence  was  picked  to 
increase  the  angle  of  attack,  the  perturbation  could  as  well  reduce  the 
maximum  angle  of  attack. 

The  wind  tunnel  test  at  AEDC  (Reference  2)  consisted  of  measuring  the 
damping  (C  +  C^)  and  overturning  (CmQ)  moments  using  a  one-degree-of 

freedom  flexure  balance.  The  wind  tunnel  model  had  slots  which  mount  the 
0.080-inch  thick  Datasheet-E  explosive  strip.  The  strips  give  various  total 
impulses  up  to  a  maximum  of  0.294  lb/sec.  The  angular  position  of  the 
model,  0,  was  calculated  from  strain  gage  measurements  taken  at  about  3500 
data  points/sec.  These  data  were  plotted  automatically  after  each  shot.  The 
predicted  angle  was  plotted  on  the  same  plot.  The  predicted  angle  is  cal¬ 
culated  from  the  measured  thrust  off  aero  data  (C  +  C  .  ,  C  ),  the 

mq  ma  ma 

measured  moment  of  inertia,  the  balance  stiffness,  and  the  statically  meas¬ 
ured  total  impulse  which  varies  less  than  5  percent  from  shot  to  shot. 

The  Schlieren  photographs  (frame  rate  7500/sec)  show  that  the  explosive 
strip  forms  a  strong  shock  wave  (Figure  28)  which  appears  to  move  out  in  a 
shperical  manner.  The  wave  slows  as  it  moves  out  and  then  drifts  back  with 
the  Mach  3  flow.  The  disturbance  appears  to  be  strong  enough  to  have  an 
appreciable  effect  on  the  missile  aerodynamics  while  it  is  acting.  The 
pretest  analytical  study  showed  that  the  time  the  disturbance  persists  is 
an  important  factor.  Using  the  frame  rate  (7500/sec),  and  counting  the 
frames  during  which  the  perturbed  flow  is  present,  gives  a  perturbation 
time  of  about  1  to  2  milliseconds.  From  Figure  27  and  the  analytic  study, 
this  is  not  a  sufficient  time  to  cause  a  noticeable  change  in  the  maximum 
angle  of  attack. 

Figure  29  shows  the  predicted  and  measured  angles  of  attack  for  one 
shot  condition.  The  plot  shows  little  difference  between  the  predicted  and 
measured  motion.  Considering  all  shots,  one  can  conclude  that  the  explosive 
blast  wave  and  associated  disturbed  flow  has  little  effect  on  the  missile 
motion  for  these  test  conditions. 

Data  were  taken  to  determine  the  effect  of  the  exposed  explosive  slots 

on  the  aerodynamics  of  the  missile.  C  and  C  +  C  .  were  measured  for  no 

ma  mq  ma 

slot  exposed,  one  slot  exposed,  and  all  slots  exposed.  Test  results  indicate 

little  change  in  C  and  C  +  C  .  for  the  different  conditions, 
ma  mq  ma 
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Maximum  Angle  of  Attack  (Degrees) 


Figure  27.  Maximum  Angle  of  Attack  versus  Perturbation  Time 
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Figure  28.  Sequenced  Pictures  Explosive  Blast  Wave 


(Degree) 


SECTION  VIII 


SUMMARY  AND  CONCLUSIONS 

The  feasibility  of  an  explosive  impulse,  laser  guided,  body  fixed  (pur¬ 
suit  guidance),  fin  stabilized  rocket  was  studied.  A  method  of  analyzing 
the  system  was  set  up  and  many  important  principles  affecting  the  accuracy 
of  the  system  were  determined. 

Two  computer  programs  were  established  to  study  the  system.  One  of  the 
computer  programs  simulated  simple  pitching  motion  of  the  rocket  and  the 
other  program  simulated  the  complete  six-degree-of-freedom  guidance  problem. 
The  pitching  motion  program  gives  the  pitch  frequency,  the  angle  of  attack, 
and  the  change  in  the  velocity  vector  when  a  thruster  is  fired  as  a  function 
of  the  thruster  characteristics  and  the  physical  and  aerodynamic  properties 
of  the  rocket.  The  6-D  guidance  program  simulates  the  entire  guidance  pro¬ 
blem.  The  program  has  variable  input  for  the  rocket,  atmospheric  conditions 
including  wind,  initial  conditions,  target  array  and  target  motion.  A  sta¬ 
tistical  summary  of  the  accuracy  on  multiple  target  arrays  is  part  of  the 
output  of  the  6-D  program. 

It  was  determined  that  to  obtain  sufficient  accuracy  for  this  system  it 
was  necessary  to  phase  the  pitch  frequency  with  the  laser  pulse  frequency. 

The  pitch  frequency  must  be  slightly  greater  than  1/2  the  laser  pulse  fre¬ 
quency.  A  pitch  frequency  below  1/2  the  laser  pulse  frequency  results  in 
resonant  buildup  of  the  angle  of  attack  and  large  miss  distances.  A  gradual 
increase  in  miss  distance  occurs  as  the  pitch  frequency  is  increased  above 
1/2  the  laser  pulse  frequency.  This  is  due  to  pitch  and  laser  pulse  fre¬ 
quency  being  out  of  phase.  The  rocket  can  be  designed  to  accommodate  pitch 
frequency  changes  resulting  from  nominal  changes  in  aircraft  launch  velocity 
without  a  large  loss  in  accuracy. 

Approximate  formulas  were  derived  to  give  the  maximum  expected  correction 
from  a  specified  rocket  configuration.  These  formulas  can  be  used  in  pre¬ 
liminary  design  to  determine  if  the  configuration  can  correct  sufficiently 
far  from  the  unguided  trajectory. 

The  relation  between  the  accuracy  of  this  system  and  the  trim  misalign¬ 
ment  was  determined.  A  simple  equation  gives  the  trim  angle  as  a  function 
of  trim  misalignment,  the  pitch  frequency,  and  the  roll  rate.  The  6-D 
guidance  program  showed  that  increasing  the  roll  rate  (to  reduce  the  trim 
angle)  from  1.2  to  1.8  times  the  pitch  frequency  improved  the  accuracy  by 
a  factor  of  10. 

Target  motion  and  wind  cause  large  trailing  miss  distances  for  this  as 
well  as  other  pursuit  type  guidance  systems.  For  the  baseline  configuration 
considered  here,  miss  distance  goes  from  about  1.8  feet  for  stationary  targets 
to  about  17.2  feet  for  target  motion  of  30  ft/sec. 


The  CEP  for  the  baseline  configuration  was  computed  for  error  sources  of 
gravity  sag,  trim  misalignment,  wind,  and  target  motion  for  a  slant  range  of 
11,000  feet.  Nominal  values  were  assumed  for  the  error  sources  and  a  partial 
operational  correction  (by  using  the  laser  designator)  was  assumed  for  gravity 
sag,  wind,  and  target  motion.  The  resulting  CEP  was  5.6  feet. 

An  AEDC  wind  tunnel  test  was  run  to  determine  if  the  explosive  impulse 
caused  a  blast  wave  or  flow  disturbance  that  adversely  affected  the  aero¬ 
dynamics  of  the  rocket.  Test  results  indicate  that  the  large  explosive 
blast  wave  and  associated  perturbed  flow  from  the  terminal  guidance  explo¬ 
sive  strip  have  little  effect  on  the  missile  oscillatory  motion.  This  is 
probably  related  to  the  short  time  the  disturbance  persists.  These  test 
results  indicate  that  the  explosive  action  time  may  be  extended,  if  desired, 
without  adverse  effect  on  the  guidance  system. 

The  initial  analytical  simulations  indicate  that  the  laser  guided  explo¬ 
sive  impulse  correction  guidance  method  is  a  feasible  approach  for  a  low 
cost  guidance  system.  Further  development  of  the  explosive  correction  method, 
additional  analytical  studies,  and  a  system  demonstration  test  is  required 
to  assure  sufficient  accuracy  of  this  system. 
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APPENDIX  A 


ANGLE  OF  ATTACK  AND  SWERVE  FOR  PLANAR  MOTION, 
ROLL  RATE,  AND  ROLL  AVERAGED  THRUST 


1.  ANGLE  OF  ATTACK  FOR  PLANAR  MOTION 


Figure  A-l.  Missile  Pitching  Motion 

The  approximate  equations  of  motion  for  the  missile  in  the  X 
plane  are 

MV  =  -D  =  -  i  pV2sC,~F 
2  d  x 


MZe  =  Fz  C0S  9  "  Fv  sin  9  +  Mg  cos  0  +  F„  cos  0 


I  0  =  M  +  F  L 
y  y  s 


Geometry  gives 


Z  =  -V  sin  (0  -  a)  ~  -  V  (0  -  a) 
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The  swerve  force  is 


F  =  QS  (-C..  ) 
z  x  v  Nor 


(A-5) 


The  moment 


M  =  QSd 

y 


(  C  |r-+c  \  +  C  “1 

l  mq  2V  ma  2V  )  ma  ; 


(A-6) 


If  small  angle  approximations  are  made  and  small  density  factors,  such  as  f, 

P 

are  neglected  when  summed  with  1,  g/d  neglected  compared  with  —  and 

Md 

the  relative  magnitude  of  the  aerodynamic  coefficients  taken  into  account, 
then  Equations  (A-l)  through  (A-7)  can  be  combined  to  give  the  following 
differential  equation  for  the  angle  of  attack  (a). 


LF 


ci  -  K.  d  -  K  a 


fC  F 
mq  s 

2K2Md 

y 


(A-7) 


where : 


M 

V 

P 

s 

d 

Q 

I 


z  = 


Missile  mass 
Missile  velocity 
Air  density 

Missile  reference  area 

Missile  reference  diameter 

2 

Dynamic  pressure  =  l/2pV 
Missile  axial  moment  of  inertia 

Missile  transverse  moment  of  inertia 

Inertial  axes 


X  ,  Z  =  Missile  axes 

0  =  Missile  angle  Figure  A-l 

0  =  Rate  of  change  of  0 

a  =  Missile  angle  of  attack 

a  =  Rate  of  change  of  a 

F  ,F  =  Forces  along  x  and  z 
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Gravitational  force 

Normal  force  coefficient  slope 


=  Damping  moment  coefficient  slope  with  respect  to 
q  =  0  here 

=  Average  side  force 

=  Position  of  application  of  Fs  (Figure  1) 

=  Damping  moment  coefficient  with  respect  to  a 


Moment  coefficient  slope  SC^^Sa 


f  =  psd/2M 


■  V/M 


K  =  j  )  C  =  pV27rd3C  /i  8 

a  2  \  d  /  ma  ma  /  y 


Fj  -  psdd2/ Iy 

K.  =  1  (c  ♦  C  .  )  -  ~T-  (c  +  C  .  \ 

a  4  d  '  mq  ma  /  16  IV  mq  ma  ) 


The  solution  to  Equation  (A-7)  is 


a  =  e 


£a  cos  aijt  +  B  sin  uijt  j  +  a^. 


a  =  e 


A  cos  a^t  -  sin  u>^tj 

r  k*  il 

+  B  ^  sin  u.t  +  u.  cos  u.t 


A,  B  =  arbitrary  constants 


Dj  =  -  K?  ^  4  ,  pitch  frequency 


CA- 15) 


“P  =  r[r-Xl  =  - 

a  L  y  2KyMd  J 


F  L  F  C  d 

s  +  s  rcq 


v  I  '  Y 
ay  2MV  C 


ma 


(A- 16) 


At  t  =  0, 


a  =  aQ  a  =  &0  then 


A  =  a  -  a 

0  p 


■  [  40  •  r  (“0  -  “p)]/“i 


(A- 17) 


(A-18) 


For  case  aQ  =  aQ  =  0,  the  values  of  a  and  a  at  thruster  burnout  time  (t) 


are 


•,s*,  [  e>T  ( • cos  v  *  s' sin  ur )  *  >] 


■  At 

a  =  a  w.  e  sin  w.  r 

T  p  1  1 


1  +1 


K. 

a 

2w, 


(A-19) 


(A-20) 


If  aij  >  >  K.,/2 


aT  ap  |  1  -  (  cos  cot  )  e*T  | 


*  Xt 

aT  ^  “pWj  e  sin  u)jT 
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CA-21) 


(A-22) 


Now  use  Equation  (A-13)  to  give  the  value  of  a  after  termination  of  the 

side  thrust.  For  these  conditions,  A  =  a  =  a  ,  B  =  (a - £-2.^ /) 

0  t  \  t  2  f  1 


a(t2)  =  6 


(a  K.a  \ 

— - £-Usin  (o.l 

<o1  2  tOj/  1 


(A-23) 


t2  =  0  at  thrust  termination 

“T  and  6^  may  be  taken  from  Equations  (A-19)  and  (A-20)  or  (A-22) 
and  (A-20) . 

2.  SWERVE  FOR  PLANAR  MOTION 

The  differential  equation  for  the  swerve  (side  deflection)  may  be  ob¬ 
tained  from  Equations  (A-l),  (A-2),  (A-4)  and  (A-5) .  Assume  small  angles, 
cos  0  =  1,  sin  0  =  0  ;  neglect  Mg  ;  substitute  (A-4)  in  (A-2).  The 

result  is 


MZ  ~  F  -  F 
e  z  x 


(-r*“)  *FS 


Divide  by  V 


PSV  /  r  A  „  \  .  psCd  Ze  Fs 

V  ~  2M  y  CNa  CdJ  “  2M~  MV 

Assume  Z^j  V  is  order  of  a  and  «  |  C. 


Z  F 

e  PSV  r  n  ,  s 

V  ~  "  2M  LNa  MV 


(A- 24) 
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substituting  from  Equations  (A-19) ,  (A-20),  and  (A-26)  and  canceling  terms 


6  =  6  (t,  -*■  large)  = 

s  t2  2 


ROLL  RATE 


(  K? 

CM  a  ]  e^T  - sin  w.t  +  2x 

d  Na  p  I  4K  0).  1 

r  a  1 


The  spin  differential  equation  is 


.  C£p  Qsd  CU°-Si  „ 

P  -  - P  -  — - ° 


p  =  Spin 

I  =  Axial  moment  of  inertia 

x 


(A- 29) 


(A- 30) 


Llp  =  Spin  damping  moment  coefficient  derivative  (constant) 
*  Spin  driving  moment  coefficient  (constant) 

Q  =  Dynamic  pressure 


Define: 


r  -  CtoQsd 

1  "  2V  I 


_  C£6Qsd 


p  —  vv 

’  C2  =  I 


(A-31) 


The  solution  of  Equation  (A-30)  is 


C2  l  V 


+  P0e 


P  -  Pa  1  -  e 


+  P06 


(A- 32) 


^2_  _  Cl6  2V 

V  P{  "  %  d 


(A-33) 


4.  ROLL  AVERAGED  THRUST 


The  rocket  will  be  rolling  while  the  thrust 
thrust  is  applied  for  a  time  period  (roll  angle 

value  of  the  thrust  in  the  mid  angle  direction. 
Figure  A- 2  it  is  desired  to  compute  the  average 
the  mid  angle  direction. 


is  being  applied.  If  the 
4  =  |t),  what  is  the  average 

(See  Figure  A-2) .  In 
value  of  the  thrust  along 


Figure  A-2.  Thrust  Averaged  about  Mid  Roll  Angle 


|  ~  constant  ,  <t>  =  |t  ,  <t>T  =  lT  >  “  lT/2 


Fg  =  Average  thrust 


k  J Fs  cos(r-  **) 


dt 


=  i  ]-  I"  sin  |  (  j  -  t  ) 

T  A  L  V  2  /JO 


F 

s 


ll 

2 


Note  if  |t  =  it  ,  (i.e.,  1/2  revolution)  then 


APPENDIX  B 


LINEARIZED  SOLUTION  FOR  THE  MOTION  OF  A  FIN-STABILIZED  MISSILE 

WITH  TRIM  MISALIGNMENT 
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Linearized  solution  for  the  motion  of  a  fin-stabilized  missile  with 
trim  misalignment: 


X  t  i$.t  X  t  i*-t  i<f> 

*=K10e  6  1  +K2Qe  2  e  +  K,  e 


(B-l) 


Where 


§  =  Angle  of  attack  of  the  missile  in  the  complex  plane 

(see  Figure  B-l) 

K10’  K20  =  Nutat*on  ancl  precession  arms 

Aj,  X2  =  Nutation  and  precession  damping  rates 

4>1 ,  <j>2  =  Nutation  and  precession  frequencies 

i  .  VT 


t  =  Time 


0 


p  -  =  Spin  rate 

=  Trim  arm. 


r-  Total  Angle 
\  of  Attack 


Figure  B-l.  Geometric  Notation  for  Pitching  and  Yawing  Motion 


If  the  magnitude  of  the  aerodynamic  coefficients,  the  mass  and  moments 
of  inertia,  and  the  relatively  low  roll  rate  of  the  missile  are  taken  into 
account  the  damping  rates  and  frequencies  of  the  missile  reduce  to  a  simple 
form.a 


1  +  4/ 

:  ~  x  p  +  X  „  a 

1,2  ~  T~  2  "  f"Ko  •  4” 

9  \r 


1,2 


-  v: 


(B-2) 

(B-3) 


For  a  simple  type  of  trim  misalignment  where  the  zero  yaw  trim  moment 

is  QsdC  6  ,  the  non-rolling  missile  will  trim  at  an  angle  6  ,  and  the 

x  mam  m 

rolling  missile  has  a  trim  arm  given  by 


For  the  case  being  considered 


This  can  be  done  since  Ka  is  the  dominating  term  in  Equation  (B-2). 


An  assumption  has  been  made  that  for  the  small  angle  of  attack  anticipated 
the  nonlinear  magnus  moment  coefficient  will  be  small.  Although  there  is  no 
magnus  wind  tunnel  data  on  the  ATR,  magnus  data  on  similar  missiles  indicate 
this  is  so. 
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1.  THRUSTER  FORCES  AND  MOMENTS 

For  the  computer  program  the  position  of  the  target  is  established  in 
body  axis  system  by  the  equations  below. 

-1 

e  =  tan  R  /  R  (r  , 

q  z/  x  (C-l) 

-1 

e  =  tan  R  /  R 

r  y/  x 

The  computer  flow  chart  logic  Figure  C-l  establishes  the  body  fixed 
quadrant  as  seen  in  Figure  C-2.  The  sequencing  of  the  thruster  may  be  done 
in  any  number  of  ways.  The  methods  shown  here  are  for  inside  out  or  out¬ 
side  in.  The  angle  r,  is  used  in  conjunction  with  the  number  of  the  thrust 
(n  =  1,  2,  3,  -  -  N)  to  determine  the  angle  x  (chi)  which  is  used  in  trig¬ 
onometric  function  determining  components  of  thruster  force  and  associated 
moments.  For  equation  logic  see  Figure  C-3.  Symbol  definitions  are  at 
the  end  of  this  appendix. 

2.  AERODYNAMIC  FORCES  AND  MOMENTS 

The  equations  for  the  aerodynamic  forces  and  moments  in  the  body  axis 
system  are  given  in  Figure  C-4.  These  are  commonly  used  equations  for  the 
more  important  forces  and  moments  for  a  fin  stabilized  missile.  The  magnus 
moment  and  trim  misalignment  terms  are  included  in  the  moment  equations. 
Figure  C-4  defines  all  the  terms  used  in  this  Appendix. 
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Figure  C-l.  Thruster  Logic  Flow  Chart 
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(n)  (n) 


y  Forces  and  Moments 


DEFINITION  OF  SYMBOLS  IN  APPENDIX  C 


C 

DZ 


The  target  off  axis  error  in  the  X-Z  plane  of  the 
body  axis  system  (Equation  1) 

The  target  off  axis  error  in  the  X-Y  plane  of  the 
body  axis  system 

The  angle  including  the  symmetrical  distributed 
thrusters  of  a  quadrant.  See  Figure  C-l 

Angular  dimensions  of  the  dead  zone 


X 


a,  b 
o,  e 


T 


(n) 


L 


(n) 


The  cone  angle  defining  the  field  of  view  from  the 
image  plane 

The  angle  defining  the  position  of  the  thruster  to 
be  fired.  Figure  C-l 

Subscripts  referring  to  option  a  or  b.  Figures 
C-l  and  C-3 

Subscripts  referring  to  n  odd  or  even.  Figures 
C-2  and  C-3 

Constant  thrust  magnitude  for  thruster  n 

Center  of  thrust  distance  forward  of  the  eg  for 
thruster  n 


F 

z 


L,  M,  N 


Q 

d 

s 


C 


za 


Component  forces  and  moments  due  to  side  thrust  in 
the  Y  and  Z  body  directions 

Component  aerodynamic  body  forces  in  the  X,  Y,  Z 
body  directions 

Component  aerodynamic  body  moments  in  XYZ  body 
directions 

Dynamic  pressure 

Reference  diameter 

Reference  area 

Axial  force  coefficient.  Computer  tabulated  and 
a  function  of  Mach  (M)  and  a 

Y  force  coefficient  slope  3C  /38.  computer 
tabulated  (M,  B)  ^ 

Z  force  coefficient  slope  3C  /3a.  Computer 
tabulated  (M,  a)  Z' 

Spin  driving  moment  computer  tabulated  (M,  a) 
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itf'UMlftitMMI  i  ii 


i 


P.  q»  r 


=  Fin  cant  constant  input 

=  Spin  damping  coefficient.  Computer  tabulated  (M,  a) 

=  X  Y  Z  body  axis  roll  rates 

=  Magnitude  of  missile  velocity 

=  pd/2V  =  non-dimensional  spin 

2  / 

=  Magnus  moment  coefficient  3  Cm  /3v36.  Computer 
tabulated  (M,  (3) 

=  Pitching  moment  coefficient  slope  3Cm/3ct.  Computer 
tabulated  (M,  a) 


=  Damping  moment  derivative  coefficient  3C 

n 

Computer  tabulated  (M,  a) 


/a  — 
/*  2V 


2 

=  Magnus  moment  coefficient  3  uny/3v3a.  Computer 
tabulated  (M,  a) 

=  Damping  moment  derivative  coefficient  3C  /3  ^  . 

m/  2  V 

Computer  tabulated  (M,  a) 

=  Yawing  moment  coefficient  slope  .  Computer 

tabulated  (M,  0) 

=  Damping  moment  derivative  coefficient  3C  /3  | . 

n  / 

Computer  tabulated  (M,  0) 

=  Damping  moment  derivative  coefficient  3C  . 

n/ 

=  Trim  misalignment  angle 
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APPENDIX  D 


EQUIVALENT  CEP  CHART 


1  1 _ I _ I _ I _ J _ I _ 1 _ L_ 

0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

a  . 

„  nun 
K  =  -  .1- 

a 

max 

Figure  D-l.  Equivalent  CEP  Chart 


98 


INITIAL  DISTRIBUTION 


HQ  USAF/RDQRM  2 
HQ  USAF/SAMI  1 
HQ  USAF/XOXFCM  1 
HQ  USAF/XOOWA  2 
AFSC/IGFG  1 
AFSC/SDWM  1 
HQ  AFSC/DLCAW  1 
AFML/DO/AMIC  1 
AFIT/LD  1 
ASD/YEM  10 
ASD/ENFEA  1 
ASD/ENAZ  1 
AFFDC/PTS  1 
TAC/DRA  1 
SAC/LGWC  1 
HQ  SAC  (NRI/STINFO  LIB)  1 
WRAMA/MMEBL  1 
CIA/CRE/ADD/PUBS  2 
AFWL/LR  1 
AUL/AUL-LSE-70-239  1 
Redstone  Sci  Info  Ctr/Doc  Sec  2 
USA  Weapons  Comd/SAPRI-LW-A  1 
AMXSY-DD  1 
AMXSY-A  1 
AMXBR-TB  1 
Frankford  Arsenal/K2400  1 
SARPA-TS  1 
USN  Weapons  Lab  1 
USN  Nav  Ord  Lab/Tech  Lib  2 
Nav  Ord  Stn/Tech  Lib  1 
Nav  Weapons  Stn/20323  1 
Nav  Sys  Cntr/Tech  Lib  1 
USN  Wea  Cntr/Code  533  2 
China  Lake/Code  40  1 
Air  Force  Weapons  Lab/Tech  Lib  1 
Nav  Air  Sys  Comd/Code-Air  5323  1 
Off  Nav  Res/Code  473  1 
NASA  STINFO  Fac/Acquisitions  Br  1 
Inst  for  Defense  Analysis/Class  Lib  1 
The  RAND  Corp/Lib-D  1 
DDC/TC  2 
USAFTFWC/TA  1 
Nav  Weapons  Lab  1 
SARWV-RDT-L  1 
Honeywell,  Inc/Aero  5  Defense  Gp  1 
Alpha  Research  Inc  1 


Ogden  ALC/MMWM  2 
AF  SPEC  COMM  CNTR/SUR  2 
DAMA/WSA  1 
SARPA-FR-S-A  1 
US  ATOMIC  ENERGY  COMM/Hqs  Lib  1 
AEDC/ARO,  Inc/Lib/DLCS  1 
AMXSY-DS  1 
AMCR-WM  1 
Nav  Weapons  Eval  Fac/WE  1 
Off  of  Chief  of  Nav  Opns(OP-982E)  1 
Nav  Research  Lab/Code  2627  1 
Cal  Inst  of  Tech/Propulsion  Div  5 
HQ  PACAF/LGWLE  4 
USAFTAWC/AY  1 
TAWC/TRADOCLO  1 
AFATL/DL  1 
AFATL/DLB  1 
AFATL/DLY  1 
AFATL/DLOU  1 
ADTC/XR  2 
AFATL/DLOSL  2 
AFATL/DLYV  1 
AFATL/DLYD  3 
AFATL/DLDL  15 
AFATL/DLDE  1 
AFATL/DLDT  1 
ADTC/WE  1 
AFATL/DLDG  6 
AFATL/DLMA  2 
AFATL/DLMT  2 
ADTC/AD  1 
Texas  Instruments  Inc  2 
Thiokol  Corp  1 
AFRPL/MKMB  2 
AAI  Corp  1 
Aerojet  Solid  Propulsion  Co  2 
NAVWPNCEN/Code  0632  1 
Arnold  Engr  Development  Cntr  1 
Hughes  Aircraft  Co/Mail  Sta  N-77  1 
GE  Armament  Sys  Dept  2 
Rockwell  Intl/-D-183  1 
AFIS/INTA  1 
BRL/EBL  1 
AGARD/NATO  1 


99 

(The  reverse  of  this  page  is  blank) 


f 


